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Microgravity Vibration Isolation Technology: 

Development to Demonstration 

Abstract 

Carlos Maurido Grodsinsky 

The low gravity environment provided by space flight has afforded the science c omm u n ity a 
unique arena for die study of fundamental and techno lo gical sciences. However, die dynamic 
environment observed on space shuttle flights and predicted for Space Station Freedom has 
complicated the analysis of prior "microgravity' experiments and prompted concern for the viability 
of proposed space experiments requiring long term, low gravity environments. Thus, isolation 
systems cypahu*. of providing ggnifieimt improvements to this random environment have been 
developed This dissertation deals with the design constraints imposed by a cce l era tion sensitive, 
■microgravity" experiment payloads in the unique environment erf apace. 

A theoretical background for the inertial feedback and feedforward isolation erf a payload was 
developed giving die basis for two experimental active inertial isolation systems developed for the 
demonstration of these advanced active isolation techniques. A prototype six degree of freedom 
di gital active isolation system was designed and developed for die ground based testing of an actively 
isolated payload in three horizontal degrees of freedom. A second functionally equivalent system was 
twih for the multi-dimensional testing of an active inertial isolation system in a red uc ed gravity 
envi ronme nt during low gravity aircraft trajectories. 

These multi-input multi -output control systems are d iscu s sed in de t ail with e s timat e s an 
pniy floor performance as well as the atntl puf«i MiK t ittf l utlira data. The 
attenuati on performance is also given for both systems demonstrating die advantages between inertial 


r 



and non-inertia] control of a payload for both the ground base envi ronment and die low gravity 
aircraft acceleration environment. 

A future goal for this area of r esearch is *r> validate itw wrimri i ap proaches «t^ 

0.01 Hz regime by demonstrating a functional active inertial fccd fo rwar d/ feedbtt* isolation sysfcm 
tforing orbital flight. A NASA IN-STEP flight e xpe rim en t has been p ropose d to a renmpii* this 
goal, and the expected selection for the IN-STEP program has been aet for July of 1993. 
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I. Introduction 

1.1 "Microgravity" Facts and Fiction 

Prom die beginning of die space program , I nma n s have b e en intrigued by die phenomenon 
nf weightlessness. Development of more accurate theories for physical and c hemica l laws have been 
possible through die study of experimental results when die effects of gravity can be sub s t ant i all y 
reduced or eliminated A permanent weightless research environment will one day be provided 
within the orbiting Space Station Freedom. However, Cor many scientific app lic ati on s, earthbound 
methods are being employed to gwnniaip. die idealized condition of weightlessness with die aim of 
improving upon current theories. 

A variety of tools are available to provide weightless conditions for materials , fluid physics, 
and basic physical phenomena research. Numerous freefall methods have been developed to provide 
short duration periods of near- weightlessness , or "microgravity". Drop tubes and towers have 
produced "microgravity'' duration's of one to ten seconds. Aircraft flying parabolic trajectories have 
generated periods of low gravity for about 15 to 20 seconds. Sounding rockets and, of course, on- 
orbit laimrh vehicles have also been used. Each technique has benefit and liability tradeoffs 
associated with it such as available lest duration versus cost per test Also, all techniques have 
specific limitati ons on die degree of ideal weightlessness which they provide. Each facility is 
typically categorized by the duration of DC or “Quasi-Steady" gravity levels achieved. However, each 
facility has vibro-acoustic responses which typically are well above a microgravity and, more 
specifically, manned orbital fac ilities have acceleration magnitudes at low frequencies well above any 
earth bound laboratory acceleration environments [1]. This section will focus on die long duration 
on-orbit facilities using measur ements of die shuttle acceleration environment and estimates of die 
future Space Station Freedom environment. 
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Interest in vibration isolation for microgravity exp e riments has incre as ed within die 
microgravity science community as the space shuttle flight program has progressed and the mail, but 
significant, levels of residual a c ce le rat i o n on die shuttle have become more widely re c o g nise d and 
documented [23]. these residual accelera ti on s result form several sources characteristic of die 
orbiting carrier and the orbital environment Very-low-frequency (dc to 10*3 Hz) ac c e lera tions due to 
drag, tidal effects, and gravity gradients contribute micrograyhy accele rat ion levels, g c lO^gp where 
g 0 = 9.81 m/sK Orbiter thruster activity can co ntrib ute 1(H to 10*2 g/g^ accelerations with 
sig nificant duratinn, but these can he predicted and controlled. The most significant and troublesome 
contribution to most experimental payloads is die moderate-frequency (10*3 to 10.0 Hz) dynamic 
spectrum of accelerations having magnitudes in the range 10*^ to 10*2 g's. This dynamic background 
is primarily due to random excitations from manned activity on die obiter as well as small thruster 
firings for orbit keeping maneuvers. However, orbiter structural dynamics and flight systems also 
contribute observable intermittent and resonant accelerations to the environment as die artnter 
interacts with its dynamic mechanical and thermal environment [4]. To categorize die disturbances 
which are present in the space shuttle and which will be present on Space Station Freedom, the 
accelerations are grouped into three frequency ranges: (1) quasi-static external disturbances. (2) low- 
frequency vibrations, and (3) medium- to high-frequency disturbances. The first category includes 
aerodynamic drag, gravity gradient effects, and photon press ur e accelerations acting directly on die 
orbiter. The second category includes excitations of the large flexible space structures doe to crew 
motion, spacecraft attitude control, and robotic anus. The diird category includes disturbances due to 
on-board equipment such as pumps and motors having a frequency range of about 10 Hz and higher. 
The range of accelerations observed on several shuttle missions or estimated for die accessible artiit is 
given in Table 1.1 [233]- 
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Source 

Acceleration, 

S/go 

Frequency, 

Hz 

Quasi-steady or constant cycle 

Aerodynamic drag 

10-7 

0 to 10-7 

Light pressure 

10-*- . 

0 to 10“ 3 

Gravity gradient 

10-7 

0 to 10-7 

Periodic 

Thruster fire (orbital) 

2X10-7 

9 

Crew motion 

2X10-7 

5 to 20 

Ku-band antenna 

1 

o 

X 

CN 

17 

Nonperiodic 

Thruster fire (attitude) 

10-« 

* 

1 

Crew pushoff 

10-* 

1 


Table LI: Acceleration Disturbances 
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The evolution of the freedom Station design has led to pvwiH»i Kmimtifig on lon g-term, 
low-gravity expe rimen ta tion in this environment. It is now obvious that most of the true 
"microgravity' experiments will require notation from this tm&MB mflli -g c an c a mtat if 
rqvoducible and useful results are to be expected. The Spree Station freedom wffl have many of the 
same excitation sources as are present on the Shuttle. However, SSF win be a wrii larger, more 
flexible structure and, hence , will be more modally dense and more easily taunted The p re s en t S pace 
Station Design has over 100 modes below a Hertz which win not only complicate die isolation of 
these low frequency disturbances but win also dirtatf. an acceler ation environment with more erejy 
at lower frequencies t han presently exhibited on the Shuttle ortnter. The low frequency regime for 
many microgravity expe rim e nt e r s appears to be die most critical a ccor di ng to die g-jitter sensitivity 
analyses which have proliferated over the past several years [6]. 

To spe ci fically address the measured acceleration environments on-board die Space stimtif- 
and expected on the Space Station freedom, a plot is presented showing areas of rr pr t rd m«gnimd<^ 
of accelerations in certain frequency bands. Figure 1.1 from a NASA SSF study reprwwi n tc the 
expected levels of excitations which wfll be present on Freedom. These ma gnitude and frequency 
depiction's of disturbances do not represent the structural response of die freedom Station or the 
Space Shuttle to expected forcing functions. However, die expected distirrhanf^g do fall above the 
"microgravity" requirement curve shown, superimposed on this plot. The measured wxykgatywi 
environment of die Space Shuttle falls below this fine in Figure 1.1, however measurements have 
shown that there is arbiter response at certain fundamental frequencies between 3 and 10 Hz [7]. The 
fundamental frequencies expected for die freedom Station are well below 1 Hz which mt expected for 
such a large flexible structure which leads to concern about the experimental ai w m mn wit 



CREH, ATTITUDE CONTROL, ETC 



Figure LI: SSF Potential Disturbances 
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12 Motivation for Active Vibration Isolation Solutions 


The evolution of die Space Station design and die documente d Space Shuttle accelerati on 
environment has led to potential limitations on long term, low gravity experimentation in these 
environments. & is now obvious that most of the true miciogravity experiments will require isolation 
from this random milli-g e nvir o nm ent if reproducible and useful results are to be expected. An 
effective isolation system must be designed with a resonance frequency below 1/V2 times die lowest 
excitation frequency of interest. Because a large part of die transient disturbances have a frequency 
range from milli-Hz to 1 Hz, it is extremely difficult to design passive isolation systems for this 
environment This difficulty is due to the absence of materials which have useful ranges of both tow- 
modulus (providing low frequency) and appropriate damping (to avoid large amplitude oscillation). 
Two stage passive isolators can decrease the frequency range, however limited damping leads to 
potentially unstable systems in die random excitation environment [8]. 

Passive isolation systems require extremely low stiffness for the isolation of small 
disturbance frequencies for typical values erf mass associated with miciogravity space experiments. In 
contrast when there are direct disturbances to a payload a small value of stiffness is not desirable 
because of space limitations to accommodate large displacements. Thus, there is a trade off, and an 
effective design would need to compensate far both direct disturbances, if present and low frequency 
base disturbances. Active systems offer significant advantages over passive systems in die orbital 
acceleration environment This is doe to the extremely small stiffn e ss e s needed to isolate against 
such low frequency base disturbances and the added capability needed to adapt to direct dis t u rbances 
for die optimal isolation of a payload. In addition, since die responses to these two excitations require 
conflicting solutions, a closed loop system is required for die inertial control of die isolated payload 
PUO]. 
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Active systems require «««"»£ of motion or position, and a feedback or feedfor war d centred 
loop, or both, to counteract mechanical excitation and to minimize motion of an isolated body. Such 
systems introduce die complexity of a high- gam control system, but offer significant advantages in 
versatility and performance [11]. To achieve a broad spe ct r um of isolatinn, both feedforward and 
feedback control loops are used r e fe ren ci ng the isolated payload to an inertial frame rather than to die 
dynamic support. 

1J Uniqueness of Active Inertial Isolation Research 

The purpose of this thesis is to complete and expand on die development of a general 
theoretical approach to die formidable problem of isolating against low frequency disturbances and 
maintaining accelerations below extremely low levels. The primary contributions of this work were 
to develop die general inxtersmnding of active inertial feedback and feedforward control of an isolated 
payload and to further this understanding by developing a six degree-of-freedom active vibration 
isolation system and experimental hardware in order to demonstrate these theoretical inertial isolation 
control techniques. This research will lead to the use of digital active control of dependent multi- 
degrees of freedom in order to inertially reference a payload needing an extremely stable platform. 
The developmental and demonstration hardware developed constitute two six degree of freedom 
testbed active isolation platforms with the unique capability of demonstrating and verifying multiple 
control approaches for die low frequency isolation problem. These multi-degree of freedom platforms 
give die control designer die availability of twelve inertial sensors and six relative sensors for the 
various states which would be utilized in a classical or modem control design. These input sensors 
are digitized and can be used in what ever control algorithm is designed dependent on die isolation 
requirements generated for a specific attached payload. The two systems developed were a six degree 
of freedom 12 bit laboratory prototype and a high fidelity 16 bit low gravity aircraft demonstration 
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system. A classical active gpring-mass-damper control and an menial feed fo r w ard nmlti-degree of 
freedom controller were demonstrated and c ompar ed to the theoretical predictions for such an 
approach. The inertial feedback compensation techniques were not demonstrated in six degrees of 
freedom because of the straightforward n ature of the controller and will be used at a later date for the 
full verification of a hybrid control approach during a p roposed NASA Office of Advanced Concepts 
and Technology (OACT) In-Space Technology Experiment 

1.4 An Overview of the Active "Micrograv i ty* 1 Payload Isolation Field 

In 1987 the Microgravity Science and Applications Division, of die Office of Space Science 
and Applications, organized Advanced Technology Development (AID) programs in order to 
develop the enabling technologies needed for die use of Space Station Freedom as a viable 
microgravity experimental platform. One of these development programs was a Vibration Isolation 
Technology (VIT), ATD. This technology development program grew hreanse of increased 
awareness that the acceleration disturbances present on die Space Transportation System (STS) 
obiter are detrimental to many microgravity experiments proposed for STS, and in die future. Space 
Station freedom (SSF). During the beginning of the VIT ATD a large gap existed between die 
necessary acceleration requirements for many experiments and die actual acceleration environment 
achieved on the Shuttle ortatex. A VIT workshop was held to bring the users together with industry 
and nnivereity technologists to establish a dialogue between the two groups to better define needs and 
requirements. The results of die technology requirements definition phase were then used to focus 
vibration isolation technology development on die critical needs of vibration sensitive microgravity 
experiments. 

An international Vibration Isolation Technology for Microgravity Science and Applications 
Symposia and Workshop was held in April of 1991 as a culmination of die technology requirements 
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definition phase of the VIT project. This Symposia/worksbop was to p r esent id evalnate the 
of the various international organizations involved in VTT and to fnmwilute plans to develop 

mutually beneficial cooperative efforts in the pursuit of VTT technology developments. This 
workshop was also used to evaluate and shape NASA's future effort s in the VIT area, specific to 
MSAD's technology involvement [12]. 

A summary of die various groups and their contributions to the field follow with their 
respective innovations and primary accomplishments. Internationally there are two groups outside of 
the U. S. who have done research in the area of "microgravity’ vibration isolation, the European 
Space Agency id NASDA. The European's began a research effort at die University of North 
Wales. D. L Jones, A. R. Owens, and R. G. Owen developed a Microgravity Isolation Mount 
(MGIM) which is a Columbus rack-based facility designed to provide active vibration isolation for 
sensitive experimental payloads. The control basis for this wok was done using die feedback of 
relative position id velocity to control a payload through a soft suspension system. This work was 
experimentally verified in only one dimension id again used no mwiwi feedback or fee dfo r w ard 
stabilization of the payload for the base cancellation of disturbances or die ability to control direct 
disturbances should they exist The NASDA pro gr a m in microgravity vibration isolation is very 
extensive and can be summarized in the general study of both passive and active devices to improve 
the microgravity environment. They have studied the use of acceleration feedback id have dene low 
gravity aircraft tests in three degrees of freedom with some limited success. [12] 

The national arena consists of basically four gro u p s who have fabricated and demonstrated 
experimental hardware for die active isolation of sensitive payloads specific to the orbital 
"microgravity’ environment. These gr o ups consisted of Honeywell Satellite Systems Operation, 
Satcon Technology, and a A. VanFlotow at MTT. Each group de s i gn e d a six degree of freedom active 
isolation system with different capabilities and limitations. The Honeywell S»trfiiti» Sy stems group 
designed id fabricated a Magnetic Isolation System (MIS) based on die control of attractive 
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electromagnets referenced to relative sensors. Hie MIS was laboratory tested in 1986, and provided 
superior six-dof is o l a ti o n, wift dose to unit transmissib ility below the break fr equenc y of about 4 H z 
and fifth-order isolation roD-oflf above [13]. Hus isolation system was die fiat of its kind. H w ew, 
it was limited in low frequency isolation bandwidth and did not incorpor at e die a dv a nta ges of 
inertially referencing the isolated payload. The Satcon Technology approach consisted of a ox-dof 
Lorentz force actnated active isolation payload itfiiring a cubic gap-error nonlinear 
controller. This e xpe rim ental package was only demonstrated in dtree-dof and does not 
inertial stabilization of die payload. The Lorentz force actuators have advantages in their bi- 
directionality and linear force to control cur rent capab ility however, they are Hmitwi in force per unit 
power as compared to attractive actuators. In addition, it is difficult to design a six degree of freedom 
system having minimal cross axis coupling using Lorentz actuators. The final group at MTT, Dr. A. 
VonFlotow et aL, developed a six degree of freedom system using piezoelectric polymer ™i«i«i as 
the actuation devices. This experimental ha rd ware was only demonstrated in three degrees of 
freedom but did include die capability of direct inertial feadh** which is advantageous as di sc ussal 
in a 1989 ASME 12th biannual Mechanical Vibration and Noise conference publication entitled “Low 
Frequency Vibration Isolation Technology for Microgravity Space Experiments" [8]. However, die 
performance of die system was degraded by the force and bandwidth limitation <: of the ■rotation 
devices [14]. 



II. Theoretical Developments in Active 
Inertial Isolation 

2.1 Single DOF Inertial Isolation Approaches 

The active isolation systems desc ri bed in this section are effective at frequencies above 0.01 
to 0.1 Hz. This constraint arises not from technology limitations, bat from practical limitations on 
the stroke needed to isolate against very low frequency disturbances. Volume constraints in die 
Shuttle and in the future Freedom Station manned environments limit die stroke of any support 
system. Far example, aerodynamic drag will act on a solar-pointing station with a frequency equal to 
that of the orbital frequency (about 90 min. per orbit). Although drag is a function of the atmospheric 
conditions during a specific mission, an average g of 10 r? will be used. Thus, die distance Freedom 
Station would travel under such an acceleration would be 2(a /tab, or 1.5 m (4.7 ft), where a> = 2sc 
/(90x60) rad/sec and a = (9.81 m/sec^)x(10 r ^ g/go), assuming null initial velocity. Thus, an isolated 
payload would be forced to follow such a large spacecraft displacement, but be active in a much 
smaller region. This active region would depend an die volume constraints of a payload in die 
Shuttle or in the Freedom Station microgravity module. 

The following one dimensional baseline cases assume die use of an attractive electromagnet 
and a Lorentz force actuator, respectively, and can be analyzed as spring-mass-damper systems. It is 
assumed that the spring and damper characteristics for die attractive electromagnet and die damping 
characteristics for the Lorentz actuator are actively controlled and translated in to w*ii»tnr response 
by a control law dependent on die response characteristics desired. Using an attractive 
electromagnetic actuator, one can produce forces in only one direction. Therefore, to achieve a push- 
pull configuration one needs to nse two opposing electromagnets acting on an armature, figure 2.1 
illustrates the two general magnetic actuator configurations: die attractive electromagnet and the 
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Lorentz farce actuator. For die attractive electromagnetic actuators, the farce pro d u ced by ode 
magnet is p roportio na l to die square of die cur rent and inversely proportional to die square of die gap. 
Figure 22 diows the magnetic circuit actuator's squared d epen d en ce on current. Became of diese 
nonlinear characteristics, a bias current linearization technique is used. Thus, the bias current is 
used to produce a nearly linear control law such that, far onaD disturbances about this current, the 
control force produced can be assumed linear. In order to control this system, one must dose a 
control loop around position and velocity feedback signals with a bias aareaL. Other nonlinearities 
due to hysteresis and saturation arise between fibre and input coil current, but are not 

si gnificant with proper care of the magnetic circuit design. 

In contrast, the Lorentz actuator can produce farces bi-directionally. The force produced by 
a Lorentz armatnr is a vector quantity equal to the cross product of cur rent and magnetic field. 
Therefore, depending on the direction of current flow in the coil, one can produce a farce in either a 
positive or negative direction. Beaux of this actuator's linear dependence on control current, 
linearization is not needed, and it is open-loop stable. The Lorentz actuator thus has advantages over 
die magnetic circuit actuator, but requires more power to produce a certain farce than does the 
magnetic circuit configuration. However, the fo rces needed to control a payload in die "wdghtkss" 
environment of space are small, and this inefficiency is not as timiting as on Earth. 

The basic concept behind diese active feedback isolation techniques is to sense position, 
velocity, and acceleration, or velocity and acceleration, and then to drive an actuator 180° out of 
phase with this signal in order to cancel a disturbance to die payload. If there is knowledge about 
certain disturbances, a feedforward loop can anticipate an excitation id react without an error signal. 
Thus, die optimal isolation system far microgravity experiments to known and sensed orth er 
di s tur banc es requires inertial isolation of fee experiment using a fe cd fo r w ard/Fecdback type 
controller. Such a controDex does not circumvent the need for relative information of the payload in 
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(b) 

(a) Attractive electromagnetic idninr. 

(b) Lorentz force actuator. 

Figure 2.1: general magnetic actuator confignratkxis. 
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Figure 12 : squared dependence of magnrtk-orcmt actuator an current 

Ob is bias current; gap hj< 1 > 2 < ^3 < W) 
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order to follow foe large motion disturbances without exceeding boundary conditions (Le^ volume 
constraints). These active isolation techniques can be impiMnented by using eifoer analog or digital 
control schemes to close die feedback and feed fo r war d control loops. 


2.1.1 Formulations for Baseline Inertial Control Cases 


The responses of a magnetic circuit isolator and a Lorentz electromagnet in one degree of 
freedom can be evaluated by their transmissibility's and effectiveness in isolating against both base 
and direct disturbances. These transmissibUities and effectiveness functions are given with a brief 
description of their formulation. (Appendix B gives a detailed analysis of foe transmissOnlity and 
effectiveness formulations and their results for a variety of feedback inertial control schemes as well 
as for feedforward transmissibility functions with bandwidth error estimates.) First, the responses, or 
transmissibUities, of foe various control approaches are generated for harmonic base excitations by 
using foe active isolation system's differential equations of morion. These equations of morion are 
formulated by using Newton's first and second laws, where foe base displacement, u, is actually a time 
function, so that n = o(t). The same is implied for a directly applied force such that, in actuality, F = 
F(t). 

In order to demonstrate foe advantages of isolating a payload by die various 
feedforward/feedback controller designs a simple one degree of freedom spring-mass-damper system, 
shown in Figure 23, is analyzed where F s is a servo force proportional to certain feedback or 
feedforward inertial control signals. This force is achieved by referencing an actnatnr to the 
appr o pri ate signal of an accelerometer, or foe integrals of this signal, either attached to die support 


structure of die mass or foe mass itself. The general equation of morion for this system is: 

mx +c(x-u)+k(x-u)+F t =0. 


( 2 . 1 ) 
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(a) Relative feedback system. 

(b) Inertial feedback and feedforward systems. 

Figure 13 : physical r e p resent a tion of active isolation systems. 
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The servo force will be defined in two distinct inertial isolation control approaches to investigate die 
isolation characteristics of these feedback and feedforward controllers. These d efin iti on s wiD be 
represented as F*ff and F^ffe where F*ff is defined as a servo force referenced to die inertial 
position and velocity of die support structure giving an i nertia l feedforward controller, nd is 
/V-finwi as the servo force referenced lo the acceleration and the inertial velocity of the isolated mass 
giving an inertial feedback controller. In order lo depict the effect of these different comrol 
approaches on die transmissibility or effectiveness of die controller to reject base dtambances die 
transmissibility function for both inertial cases will be solved. Therefore, beginnin g with the general 
equation of motion and substituting die defined servo force lepre sent a tions , 


-and. 




mx +k(x-u)+G a y b x + G^x = 0 


(23) 


(2.4) 


and, 

mx+c(x—u)+k(x—u)+G fV u+G^ f u = 0, (23) 

respectively. If one assume s monochromatic motion of die base structure, one can write the 
dis placement u as Asincot, where A is die displacement amplitude. Rearranging t e rms and using die 
following definitions, die viscous damping factors £ = c/ZmtOn, £yfj, = G v jg/2mo)j 1 , Jjyff = Gyifj/M 
ODu, a direct acceleration feedback term, G« fl^ n = G«ft. and an inertial feedforward term, Gpjfffan = 
tOnp 2 ’ 

equations (2.4) and (23) can be put into vibration notation [15] as follows: 

x(l+G+)+2^(0jk+Qfr=<nlu, ( 2 . 6 ) 


and 


I+2^co ir x+tD*x = «(co* -eo^)+i<2cD,(^-^). 


(2.7) 
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If one takes Ac Laplace transforms of x and n, equations (2.6) and (27) become: 

(** (!+<>)+ 2^Q).5 + CD* )Z(j) = C0*C/(j), (28) 

and 

(s 3 +2fyoj+a>l)X(s)=U(ij[(al 09 ) 

Taking die frequency response and calculating die magnitude of the transfer faction one actives at 
die following equations, which give die frequency r esponse transfer function or transmissibQity of die 
isolated payload to a hannocic base disttcbance far both feedback sod fecdtbiwttd nwfs res pect ively. 


X(ja» 


1 

U(J<a) 

'J 



and 


( 2 . 10 ) 


xim 

U(J(o) 


H 

r V 

“J 

'M 

T 2 - 


1 

B 

sr 

4 

f* 

f 

2^ 


2 


( 211 ) 


Equations (210) and (211) depict die feedback and feedforward inertial isolation 
approaches for attennating base excitations. In die previous formulations the s pring and dampg 
dements were assumed constant. However, if an mnhflieai is med in die payload it mner Hkdy can 
not be represented by a simple stiffness and damping factor, fa nAtiHnw if flic isolated pg yjogl is not 
passively mounted but actively controlled by an attractive magnetic suspension sys tem it has the 
advantages of a closed loop active control system where die relative spring and viscous detaaos are 
variables dependent on controller loop gains. The feedback and feedforward transfer function 
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equations depict i deal conditions where die inwtiiii a c n sor s are not b andwi dth Ihniied id have 
infinite resolution, which of course is not die case in any realistic ex per im ental system. Th erefo r e, 
die actual performance will be d ic ta ted by the noise floor and band wi dth of die control system. 

To de mon st ra te die advantages of inertial feedback and feed fo r w ard active the 

previous transfer functions will be plotted with only one control tezm in die equations Hq»rripg die 
physical characteristics of die system and its response >o die attenuatio n of base disturbances. Taking 
equation (2.10) and plotting die response for simple acceleration feedback. Figure 2.4 illustrates die 
advantage of menially increasing the dynamic mass of the isolated payload thus, forcing die transfer 
function to die left or effectively lowering die comer frequency of die system. Figure Z5 thaws die 
feedback transfer function equation (2.10) with only inertial velocity feedback These corves show 
die response of the system to increased velocity feedback drtrxminrd from die integration of an 
inertial sensor signal. The advantage of active inertial damping derived fr om an initial r efe r ence on 
the payload is the removal of the resonant response, broadening and smoothing the transition between 
the low-frequency and high-frequency regions, while reducing both die transmission and die 
response, particularly in the low-frequency range of interest Hie effect of such a system for large 
values of velocity fee dbac k gain can be understood by noting dial it is equivalent to having a passive 
damper attached between the isolated mass and a virtual inertial reference. As the inertial Humping is 
increased, the isolated mass becomes more id more tightly coupled to die (motionless) ideal inertial 
reference. In other words, the stronger the viscous term the better the attwnmtinn This type of 
response is not seen in the pure suspension case because the velocity term was determined from the 
derivative of a relative position sensor. 

To describe die inertial feedforward cancellation approach one am simply analyze the 
numerator of equation (2.11) and compare it to the typical angle degree of freedom ymg mass 
damper transfer function. The feedforward terms are simply added to cancel oat die relative stiffness 
and damping terms to theoretically attempt an infinite attenuation of base disturbances. This can not 
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Figure 2.4: inertial acceleration feedback transmissibOitj curves. 



Figure 2^: inertial damping transnussflwBty cures. 
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be achieved because of noise floor and bandwidth limitations on die feedfor w ard control bops. 
However, an increase in roll off and overall attenuation can be achieved as will be demonstrated by 
the performance of the actual hardware des ig ned. 

2J2 Six DOF Formulation for Development and Demonstration Hardware 

For die development and demonstration system har d ware designs, attractive electromagnetic 
actuators were used as the control force devices. Since the attractive magnets are open loop unstable 
and must be stabilized by an appropriate amount of relative position feedback,, one first must 
demonstrate the relative feedback control approach for a simplified single degree of freedom system 
upon which die multi-dimensional case will be formulated, 

2.2.1 Generic Actuator Control Equation 

Consider the single DOF system shown in Figure 2.6, where m and kj are the mass of the 
experimental payload and umbilical stiffness, respectively. The umbilical fliffaftss kj re pres en ts a 
flexible connection between die experimental payload and die base, which may be necessary for 
various functions like the supply of power, cooling fluid, etc.. The base acceleration u r e p resents die 
acceleration of the support structure. The control force is provided by die attractive magnetic 
ac tuator . It is well known [16] that the magnetic force f is related to die cod current I as follows: 

0 . 12 ) 


where h, Po. Af, and N are air gap, permeability of free space, magnetic pale face area, and the 
number of coil turns, re spec tiv ely. 



23 




24 


If hQ is the desired air gap at static equilibrium, die coil current Ig win be gov erned by the 
following static equilibrium equation. 

I 2 

mg+*,A = G-j, (2.13) 

where A is the static deflection in die spring. 

If x(t) is die dynamic displacement of die mass measured from die equilibrium position, die 
differential equation of motion is represen ted as. 

mx+k 1 (x-u)=f(t), (2.14) 

where 

/(0 = /-"tg-Jb 1 A, (2.15) 

and 


:\2 


f-G 


V.+i) 


( h-(x-u)y 


■f-G 


(I.+i) 2 


(h„-(x-u)) 


2 ' 


(2.16) 


In equation (2.16), i is the incremental curr ent which is small relative to Ig. Since die relationship 
between f and i is nonlinear, a linearized relationship is used from die binomial expansion of equation 
(2.16) around the current i as follows: 


/ 


I 2 

G~y+2G 

K 




11 

hi 


h 2 -i+2G-fr(x-u). 


From equations (2.13) - (2.17), 


mx+ k^(x -u) = 2 G-^-i , 

n L 


(2.17) 


(2.18) 


where 


Jl 


W=^-2G^. 


(2.19) 
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Depending on the values of kj, 1^ and G, it is possible that k^f can be negative. If one assumes 

there is no umbilical the negative from die linearized electromagnetic circuit is most 

commonly refened to as the negative magnetic field fliffaess wfaich fa sce s the cont roll e r cwrent to be 
at least proportional to a relative displacement fee d back signal in Oder to produce a closed-loop 
stable system. For toe co ntrol le r design, die control input now b e c omes 2G(I 8 /h 0 ^)L 

Since the linearization approach is tfriedy valid far a small pertmhation around the 
equilibrium position, one could directly treat /(f) as the control input; and then de t emrinr die 

corresponding coQ current using equation (2.15) and (2.16) as follows: 

(/, +0 = G~ m ( A. - (x - u))(J +mg+k l A) v7 . 020) 

Hence, for the controller one could either consider die exact relationship (2.15) and (2.16) or the 
linf-arizftd system (2,19). Appendix C gives a more detailed explanation of die linearized approach 
including a viscous term and solving far the effective stiffness and damping of an attractive magnetic 
suspension system, assuming no umbilical connection. 

2JL2 Formulation of Six DOF Equations of Motion for General Isolation System 

The following two figures and the subsequent mathematical formulations are given as 
representative of hardware systems built far die demonstration of active inertial fe e dfo r war d isolation 
techni q ues . The systems built are generally represented in the following section where the actuation 
of either system was accomplished through three actuation pod locations giving die ability to actuate 
and control the isolated payload in six degrees of freedom. A specific experimental description of the 
developed isntarinn system will be given in die experimental section of this di s s e r tati on. 
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Figure 2.7: Six Degree of Freedom System 



Figure 2£ Top View of bolatioe System 


Referring to Figures 2.7 and 2.8 the xyz coordinate axes are chosen such dial they are 
principal axes of die rigid body at the center of mass of die payload. The parameters fj, fj. and fj, 
are the control forces in the appropriate dimensions for actuator pod locations one, two, and three, 
respectively. These fo rce s provide the necessary static eqmBbrinm forces as well as the ap propri a te 
negative feedback stabilization forces to counteract die negative spring rates with die 

attractive magnetic circuits. Let m r epresent die mass of die platform and let Ig, I x , and Iy be the 
moments of inertia of die platform about the z, x, and y, axes. Twelve attractive type magnetic 
actuators provide control forces at die actuator pod locations 1,2, and 3. 



27 


If x, y, and z are the dynamic displacements cf the center cf mass of die payload in a global 
xyz directions and p, and 6 are die angular displacements about the x, y, and z axes, respectively, 
die equations of motion can be ex pres s ed an die of Newton’s and Euler's laws as follows. 

Ml, + (*«1, +*«2, = /lr +A, +/* 

1 V3 * * l » V3 

+(~ ^*1* ■*‘^2^^* 2 

J3 l « J3 - 1 
rr^ p +(rk m3Jt —+k mak -)y r = 

1 A , + OU +*.* +*.»>A», = (A + ; y A>4 

/A + OUA - W> “ M -fiA - Ai> 

ifi, +0U -*.,.)AP, - (A -A)A 


. 021) - (227) 

A 

The / ac tuato r control forces are based on a linearized control law for an attractive magnetic circuit 
as defined in die previous section. In this m»ih«ni>rir»i representation of die equations of motion for 
die six DOF case die kjn> s represent die negative spring rates associated with die linearized attractive 
magnetic circuits, where for this general representation there are no umbilical connections asa lined . 
The bias forces associated with die linearization techniques account for the residual acceleration of 
gravity in die vertical dimensions and thus, as stated previously, since equations (231) - (227) 
re p rese n t die dynamic equations cf motion, all static or equilibrium terms have been canceled 
a p pr op ria tely. 

The actuator contra fORCS ™TI be defined as frjlnws, in a vertical and horizontal dimension, 
specific to their respective actuator pod locations 1, 2, and 3. Referencing Figure 27, die Hnrarizcd 
vertical and horizontal actuator forces win be defined as follows: 
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,^,+2 C..-&-V* 

"*1» "#J» 

/ I 2 

t2s i j .0/7 ^ .• 


(228) 


The actuator constants Gj v - G 3 V , and G^h - G3h- ac tu a tor relative position displacements 

v lr * v 3r ^lr * ^3r will be defined far the general case as in the equations of motion (221) - 
(227). The actuator constants for an attractive magnetic actuator were given in section 22.1. These 
constants will be referenced to each actuator pod location 1, 2, and 3 as shown in Figures 2.7 and 2.8. 

Figure 27 is a perspective view of a prototype system developed for laboratory 
demonst rations . From this view each actuator pod location can be broken into two components 
having a coordinate relation with the global axes as follows. 
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4=-*, 

J\=Z, 

V3 v 1 „ 
«, = Y. +— X* 

*2 2 * 2 ' 

J\=Z. 

i 3 =-Y,+-X t 

•3 2 ' 2 ' 




(229) -(234) 

These two dimensional actuator coordinates will be defined as local vertical and horizontal 
components . Therefore the artnatnr constants will be defined with these local vertical and horizontal 
coordinates as subscripts. Urns, the actuator constants are defined as follows: 





G 2r 


2M&V 

G lr 


2nXV 



2M&4» 

G 2M 







where the constants Af, and N are defined in section 22.1. Referring to Figures 2.7 and 2.8 die 
possible relative displacements which can be realized by die rigid body motion of die suspended 
payload or the support structure will resub in the following definitions of relative displacements from 
die equilibrium position, defined by die nominal eqmHbrinm air gaps hpjy • h^y, and hpUi - fc pth . 
Therefore, for small relative angles, die relative displacements at die actuator locations 1, 2, and 3 for 


die vertical and horizontal directions are as follows: 
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(2-Z,)+A(6,-9.) 

<i-z,)-«e,-e.)+4(p,-p,) 
(i-z,)-^ce,-e,)-i,(p,-p.) 

-C*-x,)+4W,-+.) 

(x-x,)I+o-r,)^+4(»,-t,) 

The two systems built will be distingnishcd from each other as a laboratory prototype and 
low gravity flight test hardware, respectively. The only major difference between these two systems 
was die large one g bias operational environment daring laboratory exp e rim e n ts. This condition 
dictated a large bias current in the vertical actuators to support die weight of the payload, dms giving 
a large negative spring rate requiring a certain amount of relative position feedback, an order of 
magnitude greater than the horizontal actuators. This of course was not the case during the parabolic 
flights, since for a p p r oxim ately 10-20 seconds the isolated payload was off loaded. Now that die 
closed-loop stabilization equations for the multidimensional case have been formulated one can 
begin to theoretically represent die incremental control cur r en t i as the summation of various control 
signals representing the relative feedback and inertial feedforward control approaches. 

2.23 Control Signal Theoretical Formalation 

As depicted in section 23.1, fir die signal dimensional case, a linearized closed-loop 
attractive magnetic circuit is controlled by an incremental current i typically rep rese nt ed as being 
proportional to relative displacement and velocity, giving an effective positive wiffness and a 
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damping term for the simple one dimensional case (See Appendix Q. In line with this formulation, 
die magnetic circuits for the three riimeng ruil cmc can be analyzed by »«nniing there is minimal 
coupling due to nan-symmetric stator-armature gap displacements, and that die average farce 
produced by a specific actnatnr pair is aligned with die ce n te rli ne of die attractive magnet and 
petpendicolar to die pole face of the actuator, far small angular di s plac e me nts and translational 
motion. Therefore, the incremental control current wfll be defined in die global coordinate sy ste m 
referenced to the relative rigid body motions and die inertial parameters of die support structure. 
These globally defined control cur re n ts will then be defined in the appropriate ac tuator local 
coordinates far die control of die isolated payload. The incremental control current i wfll be defined 
in the fixed rigid body coordinates as follows: 


M 


r -C,(x, -XJ-C^x, -X t )+DjX t +D M X^ 



-c r (y, -Y.)-c>(y,<)+i>y?.+Dj t 

i. 


-C r (z, -ZJ-C^z, -Z,)+D<Z m +D m Z, 

h 


-c*<e, -ej-c^ce, -ftj+ivA +z>*e. 

H 


H^(P, - P. )-C^(J3, -PJ+ DJ.+Dfi 





The previous set of equations (237) rep rese nt the control cur re n t formulations for the 

respective body fixed xyz coordinates where the control coefficients are defined as die relative 
displacement feedback gains (amps/m), die relative velocity feedback gains 

(* m P s sec ^' “d Ibe mertial feedforward gains, (ampsAn), and 

« (amps sec/m), respectively. Tins set of equations represents a typical active (dative 


feedback control where the relative effective stiffness and damping can be digitally controlled- In 
addition the inertial feedforward terms are introduced into die control equations to attempt die 
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c ancellati on of base mechanical excitations doe to the ■»« being refe renced to its wy pt t 
environment through the relative isolation parameters, which are needed in order to have a dosed 
loop gable system as well as defining the systems cqmhbriam position. If one were to define die 
control equations for inertial fe edback as in a previous far the aae dwnwignnai active »«ni— «rm 

problem, the control currents can be represented as follows: 


fi \ 


'-c,(x, -x.)-c„v, -i.y-Ej, - e,*,' 



-C,(y, -T.)-C„(y t -Y.)-EJ r -E„ y. 



-C w (i,-Z,)-C*(i, -i,)-E,z r -E.Z, 



-c*(8,-e.)-c ir (e,-e.)-£ i e,-£i9. 

w 







Here the inertial feedback coefficients and E S fl . - . - are rirfmwt as the inertial 

acceleration and velocity feedback terms having (amps-scc^M) and (amps-sep/in) nnits, respectively. 
Using die actuator constants for die actuation locations 1, 2, and 3 one can define the control 

A 

farces, /, far these locations by defining die incremental control currents in the ap pro p ri ate local 


actuator coordinates as follows; 
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(239) 

The constants C} v ^v3v ^ lhA3h ** A* actuator constants defined by die following 
relationships formulated in section 22.1. 
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(2.40) 

Using these equations die equations of motion can be pot into matrix fonn in cider to 
calculate the theoretical transfer functions. The control currents as represented in equations (239) 
have been assumed to be in general dependent on an a p prop ri ate stable controller. For die remainder 
of this dissertation it will be assumed that die control currents defined in equations (237) are used 
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giving die relative feedback/mextial feedforward controller which is experimentally demonstrated in 
the laboratory and Leaijet hardware. 

Six differential equations of motion have been developed winch describe the motion of a 
payload isolated by attractive magnetic actuators at three planar points incorporating the center of 
mass of the system. In summary these equations are (121) through (127) where die appropriate 
constants and control forces are defined in local coordinates from equations (239). ff the sx degrees 
of freedom are truly decoupled, the system re sponse win be only due to inputs in die appropriate 
directions. This however, is not hkely and there is some coupling b et ween degrees of freedom 
because of the non-orthogonality of die attractive magnetic circuits. This non-orthogonality is very 
small, however, and will be assumed negligible. 



m. Design and Development of Active Inertial 
Isolation Systems 

The experimental ha rd wa re designed and developed Cor active low frequency aaohdon 
consisted of a ferromagnetic isolated payload actuated at three pod locations. Each actuation pod 
consisted of three or four attractive magnetic actuators arranged to produce an effective push-pull 
configuration in two dimensions giving die capability of six DOT on die ferromagnetic 

platform. Figure 3.1 shows a layout of the laboratory prototype system where die isolated payload 
was of a hexagonal design fabricated from angle iron and welded together to the dimensions shown in 
Figure 3.2. 

The relative and inertial motion of both systems, (Le. die displacement of die isolated 
payload with respect to its support e n v ironm e nt and the acc e le r ati on of other (he payload or die 
support), are measured using eddy current probes and proof-mass accelerometers, respectively. These 
signals are then run through signal conditioning c ircuits consisting of filters and gain adjustments. 
The gain adjustments are introduced to take advantage of the fell scale range of die analog-to-digital 
conversion. For die relative signals, one does not want to roll off die response at dc because die 
package is required to follow die low frequency motion of die orbiter due to rattle space c ons t rain t s 
being limited to 1 to 3 cm (394 to 12 in). This can be achieved by an integral action in die feedback 
controller. However, the accelerometer signals for the fe edback or feedforward control need to be 
rolled-off at the lower frequencies because of this same limitatinn, 

The low-pass filters are provided becase the payload, naturally being mechanically 
attenuated at higher frequencies, is limited in acceleration performance by die active control sysfcm 
noise floor. In addition to the physical control conaderatkns for the band-pass and low-pass fillers 
employed, one must consider die frequency re sponse, noise floor, and linearity of die control sensors. 
The sensor response will limit the system performance and resolution. Shaping die control 
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Figure 3*1: P rototyp e Active Vibration Isolation System Layout. 





Figure 32 : Prototype System's Ferromagnetic Isolated Payload. 
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bandwidth of die system can increase overall system p e rfo rm anc e, as defined by a broad-baml 
acceleration noise floor measurement on die controlled payload. 

The mechanical design considerations for neb systems reflect typical h ardwa r e duigp 
practices for any space bom payloads. However, there are a few considerations rnwrjvwh^g the 
structural response of the payload and support s t ructure which are important for die perfor ma nce of 
any isolation system. These concerns should be addressed In take foil advantage of pe rfo r ma nce 
capabilities and limitations of the active control system. One should design the support and payload 
structure to have their str uc t ura l modes above the controller bandwidth. If there are excited payload 
inodes within die frequency range of interest, the active relative feedback loops could potentially 
become feedforward loops driving die system unstable. In addition, the noise floor performance of 
the system will be limited by the excited modes of the structure. These factors will effect sensor 
locations where co-location of the control sensors with the artnatim device is most desirable. 
Therefore, the s truc t ura l dynamics of the s up po rt and payload are critical for observability 
considerations. 

The layout of die system is also important, especially far obiter designs, all volume 

taken by the isolation system is subtracted from the total experimental vohnne available. In designing 
a six DOF concept one mnst begin by understanding the dynamic environment in which soch a 
system is to perform. Once the isolation requirements are established, die physical limits of die 
system DOF can be defined which will then dictate the actuato/payload configuration. Once these 
isolation system requirements are established, die frill system concept is limited only by die total size, 
power, and mass constraints of the orbital carrier. 

3.1 Closed Loop Control of Active Attractive Electromagnetic Isolation 

Systems 


3.L1 Sensor Descriptions and Placements 
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In order to describe the control of a generic payload rite equations of motion for a general 
attractive electromagnetic aospension system were developed and the controller cur rents d es c rib e d 
referenced to die relative and inertial rigid body dynamics of the isolation system. These control 
states must be measured by t ran s du ce rs and used as die control inputs in order to control the 
incremental current to die actuators. In order to assess die extent to which a disturbance effects the 
isolated platform or to control the platform's motion, motion sensors are mounted an the platform or 
the support structure producing a voltage output that is proportional to the motion experienced at the 
sensor location. This section is intended to give insight into the ccnoep t s underlying die use of 
motion sensors in the stabilization and control of a platform. 

There are a wide variety of sensors available for measuring motion. Sensors which produce 
voltages proportional to linear position, velocity and acceleration are common, as are their rotational 
counterparts. A detailed description of die mechanical and electrical technology employed in these 
sensors is varied, dependent on the specific sensor, and is beyond die scope of this dissertation, but it 
is important to address die issues which affect die selection of die control sensors far this application. 
Some of the important issues are: 

(1) inertial verses Don-inertial measurement 

(2) frequency response and sensor dynamics 

(3) noise floor and parasitic sensitivities 

(4) sensor placement. 

3.1X1 Inertial Verses Non-inertia! Measurement 

Motion sensors can be grouped into two general categories: (1) inertial and (2) non-inertiaL 
An inertial motion sensor provides a measurement of die motion with respect to an inertial r e fe re n ce 
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frame, while a non-inertia] motion sensor measures motion with aspect id a noo-inenial refe rence 
frame. Examples of inertial motion season include accelerometers and rate season (rate season 
measure linear or angular velocity). Examples of non-inertial aenson mohair Linear Variable 
Differential Transformers (LVDT), capacitance probes, eddy cur re n t probes, inductive probes etc. 

To demonstra te die difference bet ween menial and non-menial season, consider two simple 
applications of an LVDT. First, in the c onfigur a t ion town in Hgure 33, mass M is sup p o rt ed by a 
spring K and viscous damper C The vertical motion u, of the base causes the motion x, of to mass 
M. The LVDT in this configuration measures the relative position of mass M, and to base, x-a 
This is a non-inertial measurement tw ~ < " wy a given voltage output corresponds id a tmkpie 
dis placem ent x with respect id a fixed reference frame attached id to base. In contrast, co n si der to 
configuration shown in Figure 3.4. The quantities M, K, C, x and u are as described previously. A 
small er mass m, supported by spring k, is mounted in a bowsing, which is in turn mounted on top of 
massM. The LVDT now measures to relative displacement of mass m, and to housing, xf-x. The 
LVDT in this configuration also provides an appro xim ate measure of to inertial a cc el eration of mass 
M assuming tot VK/m» frequencies of interest. Regardless of to motion of to base, a change in 
to inertial velocity of mass M will motion of to smaller mass m, and thus a non-zero output 

from to LVDT. The housing fixture containing to non-inertial LVDT sensor and mass m, is 
functioning as an inertial sensor for the mass M. 

The difference between these two configurations is of fundamental importance in 
wahjiiratinn ap plications. Suppose that a control system is implemented in each of to two 
configurations above, such that actuation forces F g act on to mass M, and are dependent on the 
output of to LVDT. If the purpose of the controller is to hold to LVDT output at zero as to motion 
of to base disturbs to system, then in Hgure 33 to mass M wiB follow, or "track" to motion of to 
base. This behavior is undesirable, particularly when to distiubances are seismic, and originate from 
the ba se . In Figure 3.4 the output of the LVDT is zero when the mass M is stationary with respect to 
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Figure 33: LVDT sensor measnnncnt in aon-inertial co nfi g ur ation. 



Figure 3.4: LYFT sensor measurement in inertial configuration. 



42 


inertial space (experiences no acceleration), regardless of the motion of the base. Hie bubs M is m 
this case, menially isolated from die Manic motion of die base. Tins behavior is fir more suable 
for stabilization applications. Following this tine of reasoning if the boosing was attached to die base 
structure the voltage outjxil of the LVDT would be proportional to the inertial motion of the base. In 
this case if a servo force F g were to be x t Sa xacei to this ngnal and die stiffn e ss and 
coefficients are known, motion of the base could be effectively canceled, In essence the mass M 
would seem to be in motion to an observer an die base however, in actuality, the base would be in 
motion and the mass M would be menially fixed. 

3.1.L2 Frequency Response and Sensor Dynamics of Inertial Se ns or s 

Inertial motion sensors (accelerometers and i nert ia l tale devices) have performance 
characteristics which vary with the frequency of input excitation. Mathematically, the frequency 
response of a sensor is described by its transfer function. In general die transfer function is a complex 
value dependent on frequency. The magnitude and phase of the sensor response to the input 
excitation are typically plotted separately as a function of frequency [17]. The magnmvV* of die 
frequency response is typically sufficient to describe die important features of the sensor dynamics. 
The phase information may become important in die stability of a control system which employs die 
sensor in a control loop. 

The magnitnrte function is die ratio of die amplitude of die output to the amptmuy of die 
input excitation as functions of frequency. If die sensor is intended to measure velocity, then the ratio 
of output voltage to the input velocity excitation should ideally be fiat over the frequency range of die 
sensor. Therefore, if the magnitnde of die transfer function from input to output for the kobix is 
constant over the operating range, then the sensor response is said to be flat. 
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An inertial motion sensor can only exhibit a flat response over some finite frequency range. 
Above and below the operating frequency range (generally flic flat portion of the magnitnde fimetien) 
die response is attenuated, or "rolled off”. Hgure 3.5 shows a sensor magnitnfle fanction exhibiting a 
flat response over the region from about 1.0 Hz to about 300.0 Hz. On either side of the flat region 
the response is attenuated. A h i ghp a s s filter is often nsed to mfl off the r esponse at low frequencies. 
The reason for this is two fold. In tins application the inertial sensors are used to cancel out the 
excitation effect of a base motion disturbance, however, tire practicality of tins is limited to certain 
volume constraints as was discussed earlier. There fo re, the sensor response most be toOed-off at the 
lower frequencies. In addition, electrical noise in the signal conditioning circuitry is proportional to 
the reciprocal of the frequency, and tints mcreases with decreasing frequency. This eff e ct, called 
flicker noise, degrades the low frequency performance of the sensor and therefore, an effective flat 
bandpass filter must be designed in conjunction with the a pprop ri ate inertial se n so r in order to have a 
flat response of the sensor drcuilry configuration in the control bandwidth of interest [18]. The 
highpass filter truncates the response in the degraded region so that only die reliable portion of the 
sensor response appears in die sensor output. At high frequencies, inertial masses employed in 
accelerometers cannot respond "fast enough’ dynamically to measure die motion above a certain 
frequency, hi fact, an accelerometer which operates on die principle shown in Hgure 3.4 is a second- 
order spring -mass mechanical system, and has a magnitnde function which behaves like a low-pass 
filter. In essence, the magnitnde response of die accelerometer will be flat at frequencies wdl below 
the sensor natural frequency VK/M. (Damping can extend the range up to about 0.6VK/M) 

The shape of the nmgnitnde function can be complicated, containing multiple regions of 
varying slopes. The magnitude function of the accelerometer used in this application will be studied 
later and it will be shown that this function is not flat in any frequency range. However, tins sensor is 
still quite useful. As will be demonstrated in the results section of this dissertation. See Appendix D 
for a detailed description of transfer function d etermina tion. 
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Figure 3JS: Typical sensor transfer fanctioa magnitude plot 


frequency (Hz) 
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3.LL3 Noise Floor and Parasitic Sensitivity 

As m e nti o n e d earlier, an inertial motion sensor is Smiled in die range of aeration 
frequencies diat it can be used id measure reliably. Si addition, a given sensor can only measure 
excitations which fall within a Kmitftri range of amplitudes. The mnwa c ha nge in die measured 
quantity which a sensor is capable of measuring accurately is called the sensor resolution. Amotion 
sensor can reliably measure any excitation which within this range of amplitudes bounded at <w» 
end by tbe resolution of die sensor, and bounded at die other end by the amplitude which would result 
in die maximum output range of the sensor. If die excitation acting on the sensor location is loo large 
(i.e. has too large an amplitude), then the sensor will saturate, creating an o ut p ut which is usually 
equal to the maximum output voltage of the sensor. If die excitation is loo weak (Le. has an 
amplitude smaller than the resolution of the sensor), then factors such as static electric fields and 
magnetic interference which are not related id the excitation of interest (called noise) have more 
influence on die output of die sensor than die excitation. Hus situation leads to a poor correlation 
between die excitation and die sensor o utput. In order to measure large excitations, sensor 
components are generally made stifier, and die electronics less sensitive therefore, die amplitude 
range of die sensor is very important to consider in that die sensor noise floor must be at least an 
order of magnitude less than the smallest excitation amplitude of interest and that the largest 
amplitude expected is within die amplitude range of die sensor used. 

Tbe measurement of motions having very mall amplitudes is difficult. Particularly in 
stabilization applications where dist u rbances tend to be very small, it is important to quantify 
accurately, die resolution of the sensor. Most accelerometer manufactur e rs publish values of die 
sensitivity and maximum allowable impact for die various models they produce. The sens i ti v ity is 
merely die value of the magnitude function in die flat portion of the curve [19]. When there is no fiat 
portion on the curve, sensors arc often calibrated by experimentally fitting a straight line through data 
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showing the sensor voltage output as a function of the magnitude of the input. The slope of this hne 
is also the sensitivity of the sensor, and is equivalent to die definition given pre v io u s l y. The units of 
the sensitivity are output per unit input. 

It would seem logical that one could measure as low an amplitude excitation as desired, as 
long as the sensitivity of die sensor were high enough. Unfortunately, even an accelerometer with a 
very high sensitivity win fail to provide an accurate measurem ent when the excitation becomes 
sufficiently small due to its noise p er fo r ma nce. In addition, the amplitude at which the measurement 
breaks down varies with frequency. What is needed is a plot of the amplitude at which die 
measurement provided by a sensor breaks down versus frequency. This function is c a lled the 
incoher ent power spectrum, and is obtained using signal p roces s ing techniques [20]. Ally die signal 
processing concepts required in die analysis of interest will be m e n tio ne d. 

The incoherent power spectr u m may be obtained for a given sensor using the following 
technique [20]. Two identical sensors, sensor x and sensor y, are placed in a quiet environmait. and 
are mounted rigidly to one another so dial the sensitive axes of the sensors are c oincid e nt, Records of 
sensor output versus time are obtained fra both sensors simultaneously. Autocorrelation functions 
^(x) and Ryy(x) and die crosscoraelation function Rgyfr) are calculated from die time histories of 
the output from sensors x and y. These functions are defined below, and are described in more detail 
in Appendix D. 

RjxCt) = E{x(t)x(t+ / i)} 

RyyCt) = E{y(t)y(t*t)) <3. 1) 

R xy (t) = E[x(t)y(t-Kl)} 

In these expressions, die operator E{ } represe nt s die expected value, or the arithmetic mean of die 
quantities «*clnsed in the brackets (see Appendix D). The autospectral density and cross-spectral 
density functions are cnirulntfd from the autocorrelation and e rossetar ri a tion functions using the 
Fourier transform. These functions are defined below, and are described in detail in Appendix D. 
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O b O'©) = jf «„(T)eap(-/0W)rfc 


c w (/©) = ^^?^ (t)cxp(_/Q,l)rfl 
= //^(T)«p(-7dOT)dl, 


0-2) 


where CjQtQca), GyyQco) are die autospectral densities of die autocorrelation function Cor sensors x 
and y while G^y^co) is die crosspectral density of die crosscorrelation function for die two sensors. In 
die previous expressions: 


bw = 

bandwidth of interest, oof - tt| (md/s) 

j = 

complex operator = V-l 

CD = 

frequency (radfc) 


lowest frequency of interest (iad/s) 

COf 

highest frequency of interest (rad/s). 


The functions in equations (3.1) and (32) are then used to calculate die coherence function vy®> 
as shown below. 


t2 _ fcwT 

* C„Oto)G ff (/to) 


(33) 


The coherence function is fu n d am e nt al to signal p rocess i ng. It is a cotnpiex-valiied function 
which is a measure of the agreement between the two signals, x(t)audy(t), as a function of frequency. 
Notice that if die auto-spectral density functions for die two sensor outputs, Ggx(ja>) and GjyQm) are 
equal, tea |c,0'o»f =(G„Oa>)T =|C„Cio»f. and die coherence function is equal to mrity 
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over the frequency hand. The value rf the magnitude finhww^ fimetion k flhpayw the 
legion 

os i4i si 

and a coherence ma g nit ude of rarity implies a perfect agreement between the two ngnsk x aral y. 

Suppose that a diode is imported to die configuration consi sti ng of the two sens or s m o u lted 
back-to-back. Provided that both se n so r s are functioning, they shook] both sense the nriMim and 
generate simi l ar outputs at die time of the excitation. Now sup pose that shocks similar ip die one 
described above are occurring frequently, and have pronounced frequency ramptraug between 100 
and 200 Hz. When a Fast Fourier Transform (FFT) is p e rfo rmed QO i^Mwrvtm hiflnfiw! , aw) the 
auto-spectral and cross-spectral densities are esiwiimwt from the FFT, dre auto-spectra and cross- 
spectra associated with the two sensors will both show regions of elevated energy b e twe en 100 and 
200 Hz. What is not so predictable is dial removing die shock environment does not completely 
eliminate die sensor output Small residual motion and electrical noise will cause some sensor «ign»i 

Unlike motion or an excitation which causes nwiihw outputs from both sensors, decta cal 
noise is random in nature. It is statistically unlikely that die two sensors will produce noise-induced 
voltage outputs that are die same at any given instant in time The two sensor o ut p uts due to 
excitations will tend to agree, while those induced by noise will typically disagree. The two effect s 
can be separated using the coherence function. 


The energy that both signals have in common is (subject to probability) emsnri by 

motion input, and is called the coherent power. The aAmat p o we r for both aenao r o ut p uts is 


presented below in terms of foe coherence. The coherent power for sensor 


The coherent power for sensor y is given by 

P e ,U*)=f' f G„<ja». 


(3.4) 
<3 S) 


x is given by 



49 


Notice that the coherent power is not necessarily the aame far both sensor orapots. It folkrws then 
that die incoherent power, or that energy in die si gnals by noise, is die energy that is left 

over. The incoherent power for sensors is given by 

P kI (Jo)) = (l-r l v )G ja (Jm). ' ( 3 . 6 ) 

The incoherent power for sensor y is given similarly by 

?*(/“)= (1- Y^KvOo))- (3-7) 

The areas under the incoherent power spectra are called the cumulative power s of the incoherent 
power spectra Pjcx(ja>) and P^yCeo)- 

Q*> 

•f 

* fPioU®)** OS) 

•4 

If die signal x has units of volts, then the cumulative power has units of volts^. The square root of the 
cumulative power is called the root-mean-square (rms) of the spec trum . The s pe c tr um nns of the 
signal x has units of vohs. Although there are other measures of die noise floor of a motion sensor, 
die rms of die incoherent power spectrum is a good measure of die noise floor. It is a good general 
practice to select a sensor which has a noise floor at least one order of magnitude lower Hum die 
smallest quantity to be measured in the bandwidth of interest. 

One of the problems with highly sensitive inertial sensors is that they tend to sense motions 
that they are not intended to sense. For example, a rotational velocity sensor will have «nme 
sensitivity to linear motion. This undesired response is called cross-axis, or parasitic sensitivity. The 
amount of parasitic sensitivity dial a sensor exhibits can be assessed using die coherence function. 
The coherent power is first calculated for a set of two sensors as described. Then the two sensors are 
mounted so feat' the sehsitive axes are ordbgonal to'one another. The coherent power for the sensors 
is calculated again for the sensors in this configuration. The rms of the coherent power spectra in 
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both configurations (in axis and orthogonal mounting) are calculated The ratio of the two nns values 
gives die percentage of parasitic sensitivity. A good motion sensor should exhibit very low p a r a s i t i c 
sensitivity. 


3.L.L4 Sensor Placement 


It is good practice to place die sensor as dose as possible to die pood at which a 
measurement of the motion is desired. But often, the geometry or some other constraint makes this 
difficult. In these cases it is important to be aware of the dynamic effects of placing a sensor at a 
location away from the axis of motion along which a measurement is required. 

Consider a platform as shown in Figure 3.6. The reference frame is fixed to die center of 
gravity of the platform labeled, "CG”. The positive conventions are shown for translation and 
rotation. An accelerometer is mounted at a point P, away from die center of gravity, whose position 
with respect to the center of gravity is specified by the vector r. Taking into account all six degrees 
of freedom (three translations and three rotations), die acceleration experienced at die sensor location 
P, is given by the relation [21], 

a p =a eg + (oar )+(animr ). (3.10) 

In this equation (3.10), a p is the acceleration of point P, is the acceleration of the center of 

gravity of the platform, a is the angular acceleration vector of die platform, and (D is the angular 
velocity vector of the platform. Expanded into components, die angular acceleration and angular 
velocity vectors may be described in their respective component directions as foUows: 


a=aj +a J j+a l k 
&=toj +<!>,} +®ji’ 


(3.H) 


where i,j and k are the unit vectors along the x, y, and z axes respectively. Placed at die point P, 
the accelerometer measures not only the acceleration of die center of gravity, but also the acceleration 
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Figure 3 j 6: Inertial sensor Mounted off-axis. 
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terms asso ci at ed with the rotational acceleration and the rotational velocity crossed with the vector 
position of P. The acceleration of point P can be formnlated with aD the component vector terms as 
follows: 

5 , = a* +[(®,0>, -°O r , +(<*, -(coJ+cd’K]/ 

4[^+co’)r y +(co,a),-o>,+(a,+<D/D,)r,p , (3.12) 

+{(a, +©,©>, -(a >’ +©5)r x +(©,©, -a,)r,]* 

where r x , fy, and r z , axe the components of die vector r from the eg to the accelerometer location at 
point P. If the sensitive axis of the accelerometer is oriented along die z-axis, and the accelerometer 
has low parasitic sensitivity to motion in the other orthogonal axes, then the terms in the first two 
brackets of equation (3.12) can be ignored. The terms in die third bracket most be examined carefully 
by substituting into the terms die 'worst case* values of die angular acceleration and the angular 
velocity that may occur. If die result obtained by taking die sun of all the bracketed terms is 
significantly above the noise floor of die sensor, then the tenns cannot be ignored. 

There was one type of inertial sensor used in the control systems described in this 
dissertation. This accelerometer is a proof-mass type accelerometer from Sunstrand Corporation. 
The specific accelerometer used is discussed in the next section with reference to die design issues 
highlighted in this section and some comparisons are given regarding the noise floor of die 
accelerometer used as compared to other types of inertial sensors available. 

3.L2 Sunstrand QA-2000 Accelerometer 

The Sunstrand QA-2000 transducer is an inertial, linear accelerometer which functions on 
an electromagnetic principle. Figure 3.7 provides an fllnstration of the operating principle behind 
this proof-mass accelerometer. The inertial mass is composed of a magnetized material, and dms 
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produces a magne tic fiekL An opposing magnetic field, controlled by raying die cun e nt in a coil, 
supports the inertial mass due to the repulsive force b etw een the pennancnt wwmihH ma gnetic 
fields. Hie entreat in die coil is continuously a dju s te d in a d ot ed loop, dial the repulsive 
magnetic force restores the relative «r> » i»nti reading By altering die cunent in dte coil, 

die inertial mass is held fixed with respect lo die msa case (to within tarns toferaoce). Therefore , 
die current requited to maintain <tw« condition, is at aD times p roportional to die acceleration of die 
case. The voltage drop across a resistor in die coil circak is used to determine die voltage signal, 
which is proportional to the current by Ohm's tew. 

The frequency response of this sensor type was generated as defined in section 3.1. 1.2 and is 
shown in Figure 3.8. The peaks at high frequencies were caused by die lest apparatus, and not by the 
sensor itself. The phase loss however, is caused by die sensor, but is not limiting since the controller 
bandwidth will be roUed-off. The noise floor of die QA-2000 was also measured as described in 
section 3.1.13 over a limited frequency band. The measured noise Door is presented in Figure 3.9. 

The performance goal far these vibration isolation systems in orbit is presented in Figure 
3.10. It is desirable to employ a sensor having a noise floor dial is at least an order of magnitude 
below die smallest disturbance to be rejected. Thus, if at some frequency die performance goal is 1 
micro-g, then the level of the noise floor Aould be at most , 0.1 micso-g. This requirement specifies 
a performance curve far die sensor noise. In Figure 3.10 an orbital isolation platform performance 
goal, the sensor noise goal, and the noise floors for a variety of inertial sensor is presented on die 
same plot for comparison. As shown in Hgure 3.10 die Sunstnmd QA-2000 exhibits an acceptable 
noise floor. 


3.L3 Bendy Eddy Current Displacement Probes 




og 
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Figure 3.9: Noise floor of the Snnstraad QA-2000 accelerometer. 
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Figure 3 JO: Comparison of sensor noise floors with system specification. 
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The Bendy proximity probe is a non- inertial transducer which fimrrinns on an eddy cur r en t 
principle. The proximitor, powered by a -18 Vdc to -24Vdc supply voltage, generates a radio 
frequency (if) output voltage which is applied through die extension cable and probe lead to die probe 
tip coil. The coil radiates if energy in die form of an alternating magnetic field that induces eddy 
currents in an observed conductive specimen. The eddy c ur re n ts increase in am ptirmV as the 
observed surface approaches the probe tip, thus cansing a power loss because of a bade electro-motive 
force (emf) in the probe/proximitor system and a corresponding decrease in die if voltage of the probe 
coil. As die conductive material moves closer to die probe tip, the voltage decreases proportionally to 
the displacement. The proximitor measures the peak amplitude of the if voltage to provide a 
proportional dc output voltage that is applied to associated instrumentation. If the gap remains 
constant, the output voltage remains constant; if die gap varies as a function of time, the output 
voltage changes accordingly. Therefore, die eddy current proximitor device is a dc transducer whose 
frequency response is flat limited only by die demodulation electronics. The device is not 
electromechanical in nature as die inertial transducer mentioned previously which has a specific 
bandwidth due to its physical response to the motion of its bousing. This proximity sensor measures 
the change in if amplitude due to die generated eddy currents in a conductive material so the 
frequency response is only limited by die flow of electrons in die material and therefore, is not a 
limiting factor. 

The scale factor or proportionality constant for this proximitor is generated by calibrating die 
probe voltage with a known displacement from the probe/proximitor thereby giving a voltage output 
per unit displacement This proximitor probe has a calibrated range from 50 to 550 mils and the 
sensitivity of die transducer is approximately 20 mV/miL Since die active magnetic systems 
equilibrium position and "dead band" is defined by die relative position sensors, a scale factor or 
sensitivity had to be measured for each transducer in order to normalize the state control vectors and 
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a noise floor measurement had to be made for performance predictions and system noise floor 
performance requirements. 

Unlike the QA-2000 proof-mass accelerometers, die relative probes had to be calibrated and 
normalized to each other since individual sensitivities were not available from the manufacturer. The 
probes were fastened on a fixed mounting and a conductive target was placed on a micrometer st and . 
Measurements were taken from 0 mm (probe jnst touching target) to 20 mm in increments of 0.5 mm. 
These measurements consisted of measuring the output voltage of the transducer and calibrating each 
transducer with die known sensor target position. These conversion factors for a particular probe- 
proximi tor-cable set were used to normalize die relative sensors to one another. Some different 
combinations of sets were also measured in order to determine if one part of die proximity sensor 
system was a primary source of variance in die conversion factors. The results of these calibrations 
are given in Table 3.1. As seen in Table 3.1, die largest variance between measurements for set 1 
through 6 was 0.1 V/in. From remaining measurements establishing a variance with changes in the 
probe-proximitor-cable sets it was demonstrated that a change in proximitor electronics could give a 
scale factor variance of 0.18 V/in. 


Set 

Probe 


Cahle 

Sea. QffnJ 

1 

JUNR409117 (1) 

APRR1 12981 (1) 

1 

19.57 

2 

JUNR409115 (2) 

FEBR104309 (2) 

2 

19.27 

3 

JUNR409119 (3) 

FEBR104310 (3) 

3 

19.07 

4 

JUNR409116 (4) 

APRR112983 (4) 

4 

19.09 

5 

JUNR409118(5) 

FEBR10431 1 (5) 

5 

19.31 

6 

JUNR409114 (6) 

FEBR104302 (6) 

6 

19.17 


Table 3.1 


As in the inertial transducer case, the frequency response and noise floor performance must 
be measured or known for die relative transducers in order to predict the impact on system 
performance, and to insure that the performance goals of die actively controlled platform are met. A 
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frequency response of 0 to 600,000 ipm (0 to 36 MHz) is quoted from the proximity probes 
specifications and some frequency response curves arc given in Figure 3.11 giving a number of curves 
as functions of different capacitive loads resulting from various cable lengths. The curves are plotted 
from 360 kHz in dB's showing this sensors flat response. As shown from these curves the frequency 
response is more than adequate for this application. 

The noise floor performance is important for the relative sensor because it will give a 
performance limit for the active control system on die resolution of die displacement and equilibrium 
position of the platform. This transducer's theory of operation is dictated by die amplitude 
modulation of an rf signal, due to the magnitude of die generated eddy currents p roportionality to a 
conductive specimens distance from this rf source. Hence, die noise floor performance is limited by 
the voltage source generating die earner coils rf signal. Therefore, die voltage sources stability 
dictates the amplitude stability of the probe which is independent of the conductive specimen's 
position or displacement. Since all six relative displacement probes were powered by the same -18 
Vdc supply this voltage source's power spectral density was measured. The power spectral density 
function as defined in section 3.1. 1.1, for a stationary record (stationary in time), represents die rate 
of change of the mean square value with frequency giving a measure of the noise floor amplitude as a 
function of frequency. This power spectrum shown in Figure 3.12 was plotted as the nns voltage 
verses frequency for a bandwidth from 0.1 to 100 Hz. This measurement gives die noise floor for 
each probe, in displacement units, from each probe-proximi tor-cable sensitivity factor listed in Table 
3.1. A discussion (rf this noise floors propagation into acceleration performance will be discussed in 
section 3.1.6. 

3.1.4 Signal Conditioning 


3.1.4.1 Inertial Signal Conditioning Circuits 




Fignre3.ll: Frequency response of Bentiy relative sensor. 






62 



Figure 3.12: Eddy current probes voltage source noise floor measurement. 
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In die interest of optimizing sensor range and linearity, only a relatively small portion of die 
full scale range of die accelerometer will be utilized. The closed-loop operation provides this 
capability from the QA-2000 since the torquer coil, which applies die control force to the proof mass, 
maintains the proof mass at die refe rence position while in the presence of an input force. Therefore, 
the measure of the torquer current provides a measure of the input force. Thus, a load resistor 
inserted between the current output and ground allows for the selection of the sensor's full scale range 
and provides an output voltage proportional to input force. 

This output load resistor (connection between the signal out and return pins) can range from 
a short, giving the current proportional to g, to over 10 megaohms without affecting performance of 
die transducer. Referencing the signal to ground through a smaller load resistor permits wider g 
ranges, while larger resistance loads permit sensing of snaller input forces over narrower ranges 
(higher voltage scale factors). In practice, the selection erf die load resistor can be calculated 
assuming the largest dynamic acceleration range to be sensed and multiplying this range by the 
accelerometer current scale factor, This product will be in units of milli-amperes and if one 

calls this result I, the maximum allowable load resistance Rmar (torquer resistance plus reference 
load resistance to ground) can be calculated by die following relationship: 



where V re f is the voltage source for the transducer, and Rforq is the sensors torquer coil resistance, 
0.146 m. Bom this relationship, equation (3.13) and from a desired scale factor in volts/g, K^ v ), 
one can compare the desired load resistance for the transducer circuit with the maxim um load value 
for an assumed dynamic acceleration range. 

The six accelerometers used for die determination of the six inertial degrees of freedom of 
the support structure exhibited approximately the same current scale factor, so the above calculations 
only had to be completed once. The value, &i(A) = 1-33 mA/g, along with the desired design g range 
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of ± 1 g provides the needed values far determining die load resistance for a voltage scale factor used 
in die inertial signal conditioning circuits. The following calculations support the value of a 7.5 U2 
resistor chosen for die inertial bandpass circuit design giving an approximate 10 V/g sensitiv ity 


0) I = K liA) » grange) 

I = 1.33mA/ g«lg 
/ = 1.33mA 


( 2 ) 




I 

(15-2,5)V 

L33mA 


-146Q 


=9.2 5kQ 


(3) Selected K( V) = 10 V/g 


(4) R, — K w / K I{A) 


R - 

‘ 1.33 mA/g 

R, = 7.692H2 = 7.5K2 


Therefore, selecting die load resistor as Rj = 73k£2, for all six accelerometers, a range of ± 1.0 g was 
obtained resulting in a 10 V/g voltage scale facto-. 

Unlike most quartz-based accelerometers, the QA-2000 provides a dc measure of input force. 
Thus, in the presence erf a gravity field or a large residual acceleration, the sensor responds to the 
force opposing gravity, e.g., die force resulting from the surface upon which die sensor is mounted. 



Here, this force acts on the accelerometer case and not on the proof mass, thereby causing the proof 
mass to deflect from die reference position. As described previously, die output voltage required to 
null the displacement of the proof-mass will be a direct measure of the input acceleration. Therefore, 
large bias voltages are inherent in die signals generated by this sensor and a high-pass filter is needed 
for the closed loop inertial stabilization of the isolation system as weD as the low-pass filter needed to 
roll off the accelerometer signals because of phase loss and electrical noise. Therefore, a ba n d p ass 
filter is dictated which must be designed with the closed-loop performance requirements. 

This bandpass filter was designed by cascatixmg a four pole active high and low pass filter 
with a programmable low noise operational amplifier giving a gain adjustment capability to the 
inertial signal conditioning circuit These active filters were designed as four pole Bessel fillers. 
Bessel filters give good phase characteristics, Le., constant signal delay in die passband. Which is 
very desirable where performance in the time domain is important as in this application which is for a 
near real time active control system. The amplitude response of a filter, though important is not the 
critical issue for die design of this inertial signal conditioning circuit A filter which is characterized 
by a flat amplitude response may have large phase shifts. The result is that a signal in the passband 
will suffer distortion of its waveform. Since in this application the shape of die waveform is very 
important a linear-phase filter (or constant-time-delay filter) is extremely desirable. A filter whose 
phase shift varies linearly with frequency is equivalent to a constant time delay far signals within the 
passband, i.e., the waveform is not distorted. Therefore, the Bessel filter was chosen having 
maximall y flat-time delay within its passba nd [18]. 

The specific bandpass filter design which is used for all six inertial feedforward channels is 
illustrated in Figure 3.13. The bandpass filter was designed with a “flat’ transmissibility or frequency 
response with the high pass part of the bandpass having a cat off frequency of about 0.03 Hz and the 
low pass cut off set at about 30 Hz. The 30 Hz cut off frequency for the bandpass filter was set 
because the active isolation system, being designed as a dynamically soft system, will act as a tow 
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pass filter and attenuate the higher frequencies. In addition, die payload did not need to be inertially 
controlled at higher frequencies to meet acceleration requirements. 

The inertial signal conditioning circuit has three sections. The first being a four-pole high- 
pass Bessel filter with the high-pass cut-off frequency set at 0.03 Hz, as mentioned earlier. The 
second stage having an instrumentation amplifier using a digitally controlled pro gramm a bl e gain 
Bun-Brown PGA-2000 (Programmable Gain Amplifier). The prog ram mable gain was dip switch 
selectable from 1, 10, 100, and 1000 times gain on the circuit board for each of the six accelerometer 
channels. The third stage was a four-pole low-pass Bessel filter with a cut-off frequency of 30 Hz. 
From Figure 3.13 die inertial circuit configuration also allowed for a jumper selectable by-pass 
around the active four-pole high-pass filter. This by-pass was included as an option because of the 
very slow response of the front-end high-pass active filter. If a large shock acceleration was 
introduced into the system it would take seconds to decay, causing die inertial control to react with 
excessive displacement drift The by-pass included a simple RC network (R = 1 Mil, C = 4.7 pF). 
With this by-pass selected, die front-end of die inertial circuit has die high-pass characteristic as 
follows: 


y- me 


The dc gain for this circuit assuming die instrumentation amplifier stage had a gain selection of one, 
is calculated to be 3.8. Appendix E has die detailed design and a list of the circuit component values 
and parameters. 

Following the completion of die design and assembly of die inertial circuits, a validation of 
the circuit performance and die individual circuit gains had to be measured in order to normalize the 
inertial control states to include the sensitivity factors for each inertial transducer/signal conditioning 
channel set The individual inertial circuits and a Hewlet Packard 3S66A signal analyzer were 
configured for taking data for transfer function determination. As shown in Figure 3.14, ± 15V were 








Figure 3.14: Filter transfer function determination. 
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applied to the individual filters and die filter output was connected to a channel of the analyzer. The 
tracking signal of the analyzer, which generates random noise, was connected to the filter input and 
connected to an input channel of the analyzer. The high-pass filter was not analyzed using die 
spectrum analyzer due to the long measurement time required for enough averages to determine a 
smooth transfer function. The AC coupling high-pass circuit was tested using an oscilloscope to 
determine the AC coupling ability of die circuit, as well as die instrumentation amplifiers tri mm i n g 
capability of any dc voltage offsets from internal electronics as well as the tra n s duc ers bias signal 
measurement. The transfer function was taken with ten averages in the frequency band from 0.1 to 
100 Hz. Table 3 2 gives the gain settings for a ten times gain select on the instrumentation amplifiers 
for each inertial channel. These values with the individual inertial transducers sensitivities were used 
to normalize the inertial control states for the digital control algorithm. 


Inertial Conditioning Circnit 
Vertical Location 1 Channel 1 
Vertical Location 2 Channel 2 
Vertical Location 3 Channel 3 
Horiz. Location 1 Channel 4 
Horiz. Location 2 Channel 5 
Horiz. Location 3 Channel 6 


1ft Times Cam Setting 
33.4309 
33.9174 
33.4180 
31.8346 
32.9180 
34.8583 


1 Times Setting 
335722 
3.47602 
3.45865 
327049 
3.60637 
334386 


Table 3.2 


3.1.42 Relative Signal Conditioning Circuits 


The Benlly proximitor probes used for the stable closed-loop control of the suspended mass 
have a linear range from 127 mm to 13.97 mm (0.05 to 035 in). The calibrated probes gave a scale 
factor for each probe for control channel normalization of each probe-proximi tor-cable set as 
described in section 3.1.3. These relative transducers give a voltage output proportional to 
displacement from a conductive target where die maximum voltage is equal to the proximitor supply 
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voltage source and a minimum at zero gap distance (probe touching conductive target). Since the 
platform controller is designed to suspend the mass, Le., the experimental payload, around some 
defined equilibrium condition, the relative dis pl acement output voltage needs to be referenced to a 
relative motion around this equilibrium position. 

In order to arrive at such an equilibrium state die relative controller functions as an error 
signal around die equilibrium condition defined by a relative displacement output of zero, Le., no 
deviation from a prescribed nominal equilibrium gap distance. This requirement dictates die relative 
conditioning circuits to have bucking voltage front-ends forcing the transducer outputs to have control 
voltages about zero. The relative control signals were also filtered at die higher frequencies to roll-off 
the control of the magnetic circuits. 

The relative displacement conditioning circuits consisted of two parts. The first was a 
summing circuit which acted as the bucking voltage stage to set die equilibrium position of the 
platform. The second part of the circuit was a two pole low-pass Bessel filter with a cut off frequency 
set at 35 Hz. The output was configured to be jumper selectable between the filter output or die 
summer output, (bypassing the filter). This allowed die relative control signals to have a much higher 
bandwidth, i.e., limited only by the digital loop speed, and allowing the control signals to be rolled off 
and therefore limiting the noise influencing the controllable power supplies which control the 
magnetic actuators incremental control current. 

The two pole low-pass Bessel filter was chosen because of its phase characteristics: 
a linear phase response, Le., a constant time delay in the passband, as discussed in the inertial signal 
conditioning section. Figure 3.15 is the circuit layout for the relative displacement transducers 
showing the front-end summing circuit followed by a two pole low-pass active Bessel filter. As was 
done for die inertial conditioning circuits, a validation of the circuit performance and die individual 
circuit gains were determined in order to normalize the transducer-conditioning circuit pairs to each 
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Figure 3.15: Relative sensor conditioning circuit 
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other. Transfer functions were taken for each circuit as described in section 3. 1.4.1. Table 3.3 gives 
the gain values for the six relative conditioning circuit channels. 


ttrfativ flnnititKiBmg Circuit s 
Vertical Pos. Sensor Location 1 (Sen. 1) 
Vertical Pos. Sensor Location 2 (Sen. 2) 
Vertical Pos. Sensor Location 3 (Sen. 3) 
Horizontal Pos. Sensor Location 1 (Sen. 6) 
Horizontal Pos. Sensor Location 2 (Sen. 4) 
Horizontal Pos. Sensor Location 3 (Sen. 5) 


5.14348 

5.48857 

520922 

521621 

522146 

5.18666 


Table 3.3 


3.1.5 Sensor Placement 

The number of degrees of freedom to be controlled, or control states, of the isolated platform 
were three translational and three rotational motions. These degrees of freedom were defined 
previously as the three dimensional rigid body motions of the controlled platform. As discussed 
previously «n<»rtiai and non-inertial measurements were utilized to demonstrate die active relative 
frod bar V and inertial feedforward cancellation of base induced disturbances. The accelerometer 
transducers were placed to sense the inertial rigid body motion of the support structure while die 
relative proximity sensors were placed to sense the relative modem between the controlled platform 
and the support structure. 

As dictated by the need to resolve the six rigid body motions, six equations can be written 
describing the control input channels. Therefore, six relative and inertial sensors must be poled in 
order to derive the six control states. Once a generic set of algebraic equations are solved these 
equations can be used for both the relative and inertial sensors since the placement of these 
transducers were geometrically equivalent. 
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The geometric placement of the control sensors were die same in function for both the 
prototype and demonstration hardware configurations. The exact geometrical placement for the 
relative sensors for both sets of hardware are shown in Figure 3.16 (a) and (b). The inertial 
transducer placements were located under each actuator pod, for the prototype hardware, and housed 
above each actuator pod in die Leaijet configuration. Figure 3.17 (a) and (b) show the placement of 
the inertial sensors for both hardware configurations. 

In order to derive the control degrees of freedom a generic three dimensional representation 
of the relative displacement transducers is given in Figure 3.18. This figure shows die proximity 
sensors output relation to the target displacement with respect to a global or structure fixed x, y, and z 
coordinate system. The pro ximi tors will be generically numbered therefore, die same equations will 
be applied to a generic set of inertial sensors. 



Figure 3.18 


From this pictorial representation of the transducer placements die following equations can be 
resolved for the three translational and rotational displacements of the target body from die global x. 



(a) Prototype System Configuration 
Figure 3.16: Relative sensor placement locations. 
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(b) Learjet System Configuration 
Figure 3.16: Relative sensor placement locations. 




(a) Prototype System Configuration 
Figure 3.17: Inertial sensor placement locations. 






(b) Learjet System Configuration 
Figure 3.17: Inertial sensor placement locations. 
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y, and z coordinate directions. The following equations are the global relative control states for the 
feedback/feedforward control algorithm: 


SensorS - Sensor 4 
2/* cos 30' 


(3.15) 


y= 


' Sensor A + Sensor5 
2 


■Sensor 6 


I(l+cos60')/* 


(3.16) 


z = [Sensor 1+ Sensor 2 + Sensor3]~ 


0 = [Se/tsor4 + SensorS + Sensor6]— 


(3.17) 

(3.18) 


= £j(Se/wor2 + Sensor'S) — j(S«isorl)J-^j- (3.19) 

<J> = [Sensor 3 - Sensor !]— — . (3 JO) 

2/ r sm60 

The constants l v , and 1^ represent the distance of the horizontal and vertical transducers centerline to 
a parallel line through the center of the global vertical and horizontal coordinates. These specific 
constants are not extremely important in the implementation of the control in that the controller gains 
will be representative of the constants associated with geometric positioning of the transducers as well 
as those gains fen 1 die respective proportional, derivative, and integral toms in the control equations. 

Coupling between the control degrees of freedom was assumed to be negligible far the 
theoretical development discussion of an active feedforward/feedback control algorithm. However, 
die placement of the vertical verses horizontal relative transducers will result in cross-axis stiffness 
and damp in g terms as well as the inertial transducers cross-axis sensitivity not being representative of 
input disturbances. The most important of these coupling effects is the relative transducers output 
caused by coupling between vertical and horizontal motions due to the proximitor sensors target 
cotter not being in the same plane. For instance, if their is an angular motion about an x or y axis 
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and the horizontal proximity senses' targets are not centered on this xy plane die relative transducers 
would have sensed a net translational motion in die horizontal sensors control plane defined by die 
target center points. However, this displacement was in actuality a rotation about die x or y axis and 
should not effect die horizontal control signals. This coupling can be minimized by placing the 
vertical and horizontal proximitor target center points in die same plane since this will define die 
controllers global degrees erf freedom. In practice it is difficult to realize this. However, if the 
distance between die horizontally defined center point plane and die vertically defined center point 
plane is kept at a minimum die coupling between vertical and horizontal degrees of freedom is 
minimiz ed. The prototype was an active six degree of freedom system however, since die vertical 
control outputs were significandy greater in magnitude due to die bias currents needed to support die 
weight of the platform, the horizontal and vertical control was accomplished independently. Data was 
acquired with an attempt to decouple the horizontal and vertical degrees of freedom which of course 
posses serious limitations in testing a fully functional six degree of freedom controller. Therefore, if 
there were significant rotational vertical inputs, die horizontal sensors would have sensed motion due 
to this vertical/horizontal coupling. 

In order to demonstrate the possible coupling due to this effect, assume that there is an 
angular displacement of the support structure around die x axis, such that the control algorithm 
menially cancels this disturbance and the relative modem of the platform verses the support structure 
causes the platform to go to its geometric extremes, defined by die volumetric constraints on the 
system, which in the prototype hardware was about 6.35x10*3 m (025 in.). This would be equivalent 
to an angular motion erf about 2-5x10*2 radians . This modem would be detected by die horizontal 
sensors as a net translational displacement in die y axis of approximately 1.905x10*3 m (7J5xl0* 2 
in.) because the plane defined by die horizontal sensors center targets was approximately 7.62xl0* 2 
m (3 in.) from the vertical sensor defined plane. Therefore, in order to avoid solving die six 
equations in their coupled form, the relative transducers center targets were placed in the same plane 



80 


for the Leaijet demonstration hardware, and thus, the horizontal and vertical displacement sensors 
could be assumed independent. 

32 Noise Flow Measurements and Performance Limitations 

In order id estimate the performance limitations on die closed-loop active isolation system 
one must appropriately connect die control input channels from die sensor outputs to their respec t ive 
signal conditioning circuits. The control output signal resolution was one to one with die input 
digital resolution and the noise floor of the actuator power amplifiers were found to be an order of 
magnitude lower than the digital resolution. Therefore, the noise floor performance of the system was 
dominated by the relative control inputs. The inertial feedforward cancellation approach would be 
limited by the inertial noise floor however, assuming no disturbance in die passband of die controller, 
the system acceleration performance is dominated by die equilibrium or steady state performance 
which is defined by the relative displacement control input and digital resolution. 

For a prediction of die performance limits for die relative control signals an infinite digital 
resolution, is assumed so die noise floor is defined by the proximitor/relative signal conditioning 
channel noise. Using Figure 3.12, where the proximitor voltage source stability was given as a 
function of frequency, and taking the first probe-proximi tor-cable set sensitivity scale factor from 
Table 3.1 of 19.57 V/in, one arrives at a performance limit as a function of frequency which is shown 
in Figure 3.19. The noise floor of the relative conditioning circuits was an order of magnitude below 
that of the probe-proximitor-cable noise, therefore, it does not enter into the controller performance 
prediction. The transfer function and gain of die signal conditioning circuits do affect the analog 
noise spectrum, and therefore, must be accounted for. Figure 3.20 shows die sensor channel 1 signal 
conditioning circuits transfer function illustrating die passband of this low pass filter up to the cut off 
frequency or -3 dB point at 35 Hz. The gain for this filter in die passband was measured at 1 Hz with 
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a magnitude of 5.14348. The designed gain setting was calculated to be 4 from die largest allowable 
displacement of ± 0.25 in. giving a maximum voltage output of ± 5 V from the probe-proximitor. 
Therefore, in order to utilize the ± 10 V full scale range of the digital conveners the backing pan of 
the circuits had a gain setting of 4. However, die low pass Bessel filter has some dc gain associated 
with it This dc gain value was 1.268 giving a designed value of 5.072. The discrepancy from the 
measured value of 5.14348 at 1 Hz arises from component resistance values in the circuits being no 
better than 1%. This is why each signal conditioning channel must be analyzed and the actual gain 
values measured to normalize the control input channels. From Figures 3.19 and 3.20, one can 
estimate the analog noise floor of the relative measurements to be on die order of l.OxlCT^ in. in the 
passband of the conditioning filter. The noise floor performance is a function of frequency and an 
estimate of the performance for an infini te resolution digitizer in non-dimensional acceleration units 
verses frequency can be calculated by taking Figure 3.19 and dividing it by a gain of 5.14, taking the 
second derivative in the frequency domain of this curve, and dividing by 386.4 io/s^, thus non- 
dimensionalizing the curve. This non-dimensional noise floor performance is given in Figure 3.21. 
This prediction was estimated assuming that the digitizer had an infinite resolution. Of coarse this is 
not the case and, in addition to die digital resolution or "dead band" of the digital controller (which 
will effect the system depending on die proportional, integral, and derivative (PID) terms in die 
control algorithms and the digital loop speed of the control program), die analog-to-digital (AD) and 
digital-to-analog (DA) also have some noise floor performance associated with them. Since these 
boards are tied to the same power bus, the DA card was given a command for a zero output and the 
DA was connected to a signal analyzer measuring die voltage output noise floor performance of the 
converter. This noise floor tracked that of the proximitor voltage source, therefore, the noise floor 
predicted in g's from the proximitor estimate will be representative. 

The "dead hand" of a control system typically refers to die digital resolution of die control 
system. One can calculate from a unit scale factor, for some digitized transducer, the resolution of the 
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Figure 3.19: Proximitor noise floor due to voltage source stability. 
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digitized signal depending on the converter and die full scale maximum range of the convened 
signal. The relative displacement control is realized by six displacement transducers where these 
sensors are electronically conditioned, as discussed earlier, and digitized by a 12 bit converter and a 
16 bit converter for the laboratory and Leaijet hardware, respectively. Therefore, far the Learjet 
hardware, since the scale factor for channel 1 was computed at 19.57 V/in and the conditioning 
circuit had a gain of 5.14 in the passband, a nominal sensitivity of 100.59 V/m is used to determine 
the controller with a full scale range of ± 10 V. Tims, the unit displacement resolution or controller 
"dead band" is calculated as follows: 

^ f Full Scale Range 1 (3 .2L) 

^2 16 bit resolution J Transducer Sensitivity 

where for this case the system "dead band" was 3.034X10 -6 in, while for the 12 bit laboratory system 
die horizontal degrees of freedom had a system sensitivity of 1 .526x10'^ in. 

3.3 Actuator Design and Current Control 


The attractive magnetic actuator force is controlled by the introduction of a fluctuating 
current through the magnetic circuits conductive coil, based on die platforms relative position, 
velocity, and the inertial states of the support structure. This control current must be commanded by 
the digital controller output signals produced by the digital-to-analog converters. The control voltage 
si gnal controls the current to die magnetic circuits thus, producing control forces on the 
ferromagnetic platform. 

The maximum current requirement and voltage to current sensitivity most be specified based 
on the operational range erf the actuator magnetic circuits. These current ranges are based cm the 
operational environment for the isolation system. Since die high fidelity 16 bit Leaijet demonstration 
hardware was to be tested in an off-loaded condition, die dynamic environment of the low gravity 
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maneuver dictated the control current range needed far the current power amplifiers. The maximum 
current output per actuator and the current amplifier control sensitivity were also based on the 
actuator magnetic circuit design. 

The actuator designs for the laboratory and Learjet hardware were different due to then- 
different operating environments. The laboratory system had to contend with a large 1 g bias field in 
the vertical direction, and therefore, die vertical control magnets were sized accordingly where as the 
horizontal actuators did not su pport a large bias force. These actuators were sized to control the 
horizontal degrees of freedom in die microacceleradon range. Hus section will be dedicated to the 
design and control of the Leaijet demonstration system and will be kept as a general discussion. 
Therefore, the approach can be utilized for any attractive magnetic system based on die systems 
operational environment 

Volume constraints will dictate a certain design envelope for an attractive magnetic actuator. 
However, the most critical design specification for any magnetic actuation device is the peak force 
requirement for a specific situation or operating environment Specific to die Leaijet a dynamic and 
non-stationary loading occurs during a low gravity maneuver imposing strict limitations an system 
requirements which are not limi ted on Earth or orbi tally based platforms. The mane uver is 
characterized by three phases, an entry phase, low gravity phase, and an exit phase. These three 
phases can be characterized primarily by the pitch axis accelerations imposed on any equipment 
attached to die aircraft structure. The first phase of die trajectory is defined by a high g maneuver 
pitchin g the nose of die aircraft to a desired 45° nose high condition prior to die push-over, which 
be gins the low gravity portion of die trajectory. The nominal loading of the aircraft and its contents 
at this point is typically from 2 to 3 g's. This condition lasts a few seconds where die nose is then 
pitched forward and die free-fall begins with an initial negative acceleration followed by a 2 to 4 
second oscillatory condition while die pilots are stabilizing die control point of die trajectory, "zero 
g", after which the 10 to 15 second low gravity stabilized portion of the trajectory follows. This 
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portion of die trajectory is the so called operating environment of die platform. However, since this 
condition does not simply present itself and has a prior and future history the operational 
environment for die whole trajectory must be taken into account [22,23]. 

The maximum force of die pitch or vertical actuators was dictated by the 3 g maximum jnll 
up condition. The volumetric constraints on die actuators forced the designs to be constrained by a 
ma gne tic pole face area of (0.75 in) 2 and (0.5 in) 2 , for the vertical and horizontal magne ts 
respectively. By determining die maximum force needed, a maximum bias current for each vertical 
magnet can be calculated from the following magnetic force equation (3.12). (See Appendix B for a 
detailed discussion of magnetic circuit equations.) 


F = M^ 
“ hi 


(3.22) 


Using this equation, the magnetic cores pole face area A, and die number of coil turns N available 
around the magnarir. circuit, the maximum bias current was calculated. From the dim e ns ion of the 
vertical act uato rs a packing area of 0.352 in 2 was calculated where the conductor used for this 
application had a crossectkmal area of 0347x10" ^ in 2 . Therefore, die total number of coil turns 
possible, with a 100% packing efficiency, would be die ratio of these two areas. However, since a 
100% parting efficiency is not possible the number of maximum coil turns was 750 per magnetic 
pole giving a respectable 74% packing efficiency from a jumble coil winding process. Therefore, the 
maximum coil current for one actuator was calculated with a nominal magnetic gap of 0.25 mches 
foam equation (332) as ig^ = 334 amps. The maximum force per actuator was calculated from die 
total platform and support electronics mass of 13.61 kg. 

Another secondary issue concerning the maximum residual control current which will effect 
die overall design of the magnetic actuator coil design and its interface to die overall system 
configuration is the hearin g of the conductor coil. In a space based application this residual current 
will be dictated by die bias current linearization and will not be a limiting factor. However, the 



laboratory based system needs to account for tbe heating of die vertical coils doe to the large bias 
c ur rents needed to support tbe weight of the system as well as the maximum force needed to control 
the dynamic and steady state operational conditions of the platform. The criteria for the hea t ing of 
the coil is to insure die coils can reject enough heat so die temperature does not exceed die maximum 
allowable for die conductor msnlarinn For the Leaijet application the maximum current condition 
would only last a few seconds therefore, even though this maximum current would overheat the coils 
the temperature would not exceed specified limits in such a short time. 

The power operational amplifier used to command the actuator currents was an APEX 
PA10. This power amplifier is a high output current operational amplifier designed to drive resistive, 
inductive and capacitive loads. The features of this power amplifier gave a wide supply range of ± 10 
V to ± SO V with a peak current output of ± 5 amps. Figure 3.22 shows the APEX PA10 power 
operational amplifier control circuit which was designed to give a full scale range of 4 amps for a 10 
volt input range giving the current to voltage sensitivity of 0.4 amps/Volt. The dc supply voltage for 
the APEX power operational amplifiers was given by two Power-One Model SPM3 ± 24 V @ 5 amp 
supplies wired in parallel giving a maximum current draw of 10 amps at ±25 volts. 

As was done for the signal conditioning circuits and control sensors, the control power 
amplifiers and the magnetic actuators needed to be tested in order to normalize each control output 
channel for discrepancies in power amplifier sensitivities and actuator-current-to-force sensitivities 
between c hannels to minimiz e coupling due to non-symmetric actuation of the platform. The actuator 
power amplifiers were individually tested to determine each channel's current to voltage sensitivity. 
The individual magnetic circuits were tested by connecting each magnetic pole coil in series and in 
parallel, foe configuration flown on foe Leaijet was with parallel wired coils, a current source was 
thm provided to the magnetic cofl circuits and a Hall effect probe was placed in the center of the pole 
face with foe nominal 0.25 inch gap between foe magnetic pale face and the ferromagnetic platform. 
The sensitivity of the force producing magnetic field in each magnetic circuit was measured verses 



tbe magnetic circuit's current source. The sensitivity of each magnetic actuator was used as die 
normalization for die control force magnets. Table 3.4 gives die power amplifier circuits sensitivity 
and die attractive magnetic actuators sensitivity for die three actuation pod locations for the vertical 
and horizontal magnetic pairs. 


Magnetic Circnit Location 
Location 1 VI 
Location 1 V2 
Location 1 HI 
Location 1 H2 
Location 2 V3 
Location 2 V4 
Location 2 H3 
Location 2 H4 
Location 3 V5 
Location 3 V6 
Location 3 H5 
Location 3 H6 


Actuator Scnatiritr (GJA1 
Series; 1030 Parallel: 500 
Series: 1060 Parallel: 529 
Series: 903 Parallel: 443 
Series: 917 Parallel: 459 
Series: 852 Parallel: 422 
Series: 1000 Parallel: 500 
Series: 944 Parallel: 458 
Series: 879 Parallel: 432 
Series: 943 Parallel: 462 
Series: 1020 Parallel: 494 
Series: 952 Parallel: 467 
Series: 903 Parallel: 447 


Power Amplifier Sens, f A /VI 
0.4326 
0.4486 
0.2874 
03382 
0.4375 
0.4345 
03249 
03319 
0.4900 
0.4631 
03261 
03317 


Table 3.4 
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Figure 332 : Circuit for AFFIX PA12 power operational amplifier. 






IV. Experimental Results 


4.1 Laboratory Experimental Results 

The development and proof-of-concept or feasibility of active inertial isolation was 
accomplished using an active m a gne t ica lly suspended platform. The isolation goal of a generic 
payload was to develop die capability for die attenuation of sub-Hertz frequency disturbances to 
microgravity acceleration levels. The ground based testing of six degree of freedom systems, in all 
dimensions, poses severe limitations in bandwidth and acceleration magnitude limits. The large 1 g 
bias acceleration in the vertical direction limited die prototype hardware testing to the horizontal 
plane. However, the prototype system was actively controlled in all dimensions where the vertical 
control signals were acting primarily as an off-loading device for die performance testing of the 
horizontal degrees of freedom. The first section in this chapter will describe die experimental setup 
used to acquire the laboratory performance data for the horizontal three degrees of freedom, two 
translations and one rotation about a vertical dimension. 

4.1.1 Laboratory Experimental Environment 

The laboratory hardware was structurally mounted to a large aluminum plate which was 
bolted to the comer of a pneumatic, isolation table. This isolation table was simply used as a support 
device and the pneumatic legs were not pressurized during the testing of the suspended platform 
Initially the pneumatic table was pressurized in order to actuate die support structure of die active 
platform with a minimum amount of force. However, the table would undergo large vertical rocking 
motions at 0.6 Hz because the horizontal degrees of freedom would couple to the vertical motions of 
the pneumatic table. Hence, the legs could not be pressurized during die experiments. However, die 
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voice coil shaker was s ufficientl y powerful to introduce enough accelerations into the support system 
for calculating the desired transfer functions in die horizontal degrees of freedom. The pneumatic 
table with die prototype active isolation hardware is shown in Figure 4.1. As depicted in the 
photograph the magnetically levitated platform was attached to an aluminum interface plate which 
was fastened to the pneumatic table. Not shown in the photograph, a voice coil was a ttached to the 
pneumatic table actuating against a laboratory walL This actnator was powered by a controllable 
current source where the control input was introduced by die spectrum analyzer. A biaxial 
accelerometer head was placed on the table and the horizontal axes of the triaxial h ead were 

aligned with the x and y axes of the levitated platform on which a second triaxial accelerometer was 
mounted. These sensors were used to give the noise floor performance of the isolated platform as well 
as the transfer function performance of the forced support structure in the x and y axes. The full six 
degree of freedom performance of such isolation systems will be covered under die demonstration 
hardware results since the full six degree of freedom control was representative of an off-loaded 
system, although for a limited amount of time. 

The laboratory experimental environment was dominated by the bias acceleration in the 
vertical direction. However, the prototype system could be tested in the sub-Hz frequency d omain to 
the noise floor regime of the sensitive accelerometers by exciting only die horizontal plane and testing 
the control functions of the horizontal degrees of freedom independently of the vertical degrees of 
freedom. The nan-forced horizontal accelerations of the isolation system support structure gave the 
tt-firing resolution of the laboratory environment where die disturbances to tins support were structure 
bom due to its attachment to the building. It is assumed that die direct disturbances to the pneumatic 
table were minimal. Figures 42 and 4.3 show frequency spectrums of an x and y axis acceleration 
time history for the pneumatic table where their directions where defined by the platform fixed 
coordinates. Referring to Figure 3.8, die laboratory acceleration environment is in die noise of the 
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Figure 4.1: Laboratory view of active isolation hardware. 




tuiuuOH 

MY6 16/f/fi 


94 



Figure 4.2: Acceleration spectnuns of table x, and platform x and y directions. 
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QA-2000 accelerometer signal. From Figures 42 and 4.3, a pneumatic table support structure mode 
is evident at 3 Hz where multiples of 3 Hz are also present in both die x and y sensors, respectively. 

4.1 2 Laboratory Data Collection and Transfer Function Measurements 

The acceleration noise floor performance of the laboratory system is formed by measuring 
the accelerations erf die levitated system in die horizontal dimensions. The acceleration data was 
taken with Sunstrand QA-2000 proof mass linear accelerometers housed in a triatrial sensor head 
with two of the orthogonal sensors aligned with the platform defined x and y axes, respectively. 
These accelerometer outputs were digitized using an HP analyzer where the power spectral density 
functions were calculated. These power spectral densities of the levitated platform's horizontal 
accelerations gave the acceleration noise floor performance of die system as a function of frequency 
while the support structure was not being disturbed. These measurements were taken in an open 
laboratory so their was no attempt at shielding the experimental platform from direct disturbances. 
Therefore, these measurements assumed that the accelerations introduced to die payload were 
Hrwninafpttl by the control current noise in die electromagnets, due to die total control system noise. 
The many contributors to this overall noise in die control magnets were discussed in the previous 
chapter, and will be compared to the measured values in the experimental results discussion section. 

Figures 4.2 and 4.3, showing the acceleration spectrum erf the table x and y axes, also show 
the corres ponding accelerations of die x and y sensors attached to die isolated payload. Far 
comparison the microgravity acceleration requirement for die space station is superimposed on these 
plots, in the bandwidth from 0.1 to 10 Hz. As shown in these plots the support structure modes are 
not Higtiirhing die isolated payload and die system acceleration noise floor performance is simply 
dictated by the closed loop performance of die magnetic actuators feedback loops. Again the direct 
disturbances are asgimeri to be negligible and die noise floor is dictated by the control resolution of 
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die total system. These power spectral density functions were calculated with 401 frequency lines of 
resolution for die frequency band of 0 to 12.5 Hz. Therefore, die spectral resolution for the 
bandwidth shown from 0.1 to 10 Hz was 0.03 Hz. Six averages were used for these plots with a non- 
overlapped Hanning window. The control was strictly relative feedback since the support s tru c ture 
was not being disturbed and therefore these plots represent the system equilibrium or steady state 
noise floor performance for the 12 bit laboratory hardware configuration. The relative position 
feedback for the x, y and ^ degrees of freedom were set to give a damped natural frequency of 0.6 Hz. 

The transfer function determination of the laboratory hardware was accomplished using both 
the support structure and platform mounted accelerometers. The support structure was actuated by a 
Lorentz actuator attached between a comm' of the pneumatic table and a laboratory wall. This 
actuator was controlled by taking an input voltage signal generated by the HP analyzer and driving a 
bi-directional current controlled source provided by a laboratory controllable power supply. This 
current drove the voice coil giving an excitation force proportional to foe current controlled signal. 
Since the support structure for the levitated platform was actuated by one servo-force mechanism, the 
actuator was judiciously positioned to introduce a torque into die support table in the horizontal 
plane. This input force ran<yd the pneumatic table to pivot about its support legs since the actuator 
was positioned with a moment arm to a leg attachment point This multi-dimensional excitation at 
the isolation systems attachment points gave an approximately equal acceleration spectrum in the x 
and y axis aligned accelerometers attached to the support The angular accelerations introduced by 
die torquing of die pneumatic table were not measured since only two triaxial acceleration heads were 
available, one for the isolated platform and the other for the support structure, an angular acceleration 
measurement was not possible for die laboratory hardware. 

The spectrum analyzer assumes that the transfer function calculations or frequency response 
functions are referenced to a single channel of die analyzer. In order that multiple tests need not be 
run during die proof of concept phase of development for each control parameter settings, die support 
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structures acceleration spectnxms were checked to see if the acceleration input power spectral 
densities were truly similar in their spectral content. The two table-fixed accelerometers were 
connected to the spectrum analyzer to calculate the input power spectral densities commanded by the 
excitation source. Figures 4.4 and 4 JS show the input power spectral density functions for the 
corresponding table-fixed x and y axis accelerometers. Included for comparison on Figures 4.4 and 
4.5 the x and y platform accelerometers are also plotted and the space station requirement is 
superimposed for reference. The x and y axis support accelerometers demonstrate very similar 
spectral response. These commanded excitations were introduced by controlling the servo-force 
mechanism with a white noise disturbance signal. However, since the white noise introduces a 
voltage signal across all frequencies the current source could not drive the voice coil actuator to excite 
the support structure as strongly as if a single frequency excitation was introduced as the drive signal. 
The low frequency range was of most interest and a large force was required to drive the large mass 
of the table in the testing bandwidth from 0.1 to 10 Hz. Figure 4.6 gives the input power spectrum of 
the table -fixed accelerometers for a swept sine wave excitation input with die superimposed space 
station requirement as a reference. As shown the table accelerations are of larger magnitude than in 
the white noise excitation disturbance. 

Transfer functions give a measure of the attenuation and amplification profile of an isolation 
system. These transfer functions and their associated input and output power spectral densities are 
given in the 0.1 to 10 Hz bandwidth. This bandwidth was chosen because the goal was the sub-Hertz 
isolation of a package by active inertial means. A space bom payload will most likely have a cut off 
frequency set at about 0.01 Hz. However, performing a multi-dimensional test on a ground based 
system in this frequency range would be cost prohibitive. Their have been single-axis active tests in 
the 0.01 Hz frequency range however, die active system tested was constrained to this same axis. In 
addition, these tests did not demonstrate inertial isolation, but a simple lead-lag compensator 
approach was tested. The strokes required for a multi-dimensional test would dictate a very 
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complicated and costly test apparatus. However, testing in the 0.1 to 10 Hz frequency range for 
ground based systems is a valid demonstration for active tow frequency inertial isolation concepts. 
The approaches and active control techniques are readily extendible to the lower frequencies and 
should pose no limitations on controller designs. A second levitated system with die ap propr ia te 
noise levels and strokes could be built to test a ground based system at the lower frequencies in the 
horizontal plane, but would be unnecessary and increase costs. 

To demonstrate the advantage of active inertial isolation the transfer functions were 
calculated for an inertially referenced and a relative feedback controller. These transfer functions 
were calculated using the table x axis accelerometer as the spectrum analyzers input power reference. 
Again, referring to Figures 4.4, 4.5, and 4.6 the support structure bad very similar spectral densities 
in both die x and y dimensions however, the x axis consistently had greater acceleration magnitudes 
and, therefore, was used as the input power reference for the transfer function calculations in order to 
give conservative results. Figure 4.7 gives the transfer function calculations for die x and y axis 
platform accelerometers referenced to the input power spectrum from the support structure. These 
four curves demonstrate die advantages of inertially referencing the payload through the feedforward 
cancellation of the support structure excitations as opposed to a simple relative feedback isolation 
control approach. The natural frequency for both sets of curves was set at 0.65 Hz. As shown, the 
inertial controller is attempting to cancel die coupling of die levitated mass through the spring and 
viscous damper connections, due to the active relative position and velocity feedback of the closed 
loop relative control. The relative controller demonstrates a nominal 50 dB/decade roll off while the 
mertial controller is rolling off at approximately 100 dB/decade. The theoretical prediction for the 
inertial feedforward approach for an infinite resolution controller would show the attenuation 
continuing however, the 12 bit system reaches its resolution or deadband limits and the transfer 
function goes to (he limiting case as demonstrated by the inertial curves in Figure 4.7. In addition, 
the phase and linearity of the accelerometers and relative proximitors will effect the roll off and 
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Figure 4.6: Swept sine wave acceleration spectrum for table fixed x, and y directions. 
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solid line: plate y/ table x 
dashed fine: plate x/table x 

Figure 4.7: Transfer function calculations for relative and inertial control 
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attenuation profiles. This can be shown by plotting the input and output power spectrums of the 
transfer functions and analyzing the physical parameters of the controller from the acceleration 
spectrums. The results discussion will give a comparison of die predicted deadband and noise floor 
limits with the measured experimental results. 

4.2 Lear jet Demonstration Hardware Experimental Results 


The six degree of freedom demonstration hardware was flown through low gravity Keplarien 
trajectories to acquire performance data in an off-loaded environment Although the low gravity 
environment is limited in time and the non-stationary aspects of die maneuver cause limitations in 
bandwidth and system control parameter testing, this environment allowed die testing of the full six 
degrees of freedom with comparable control and equilibrium stales for both vertical and horizontal 
motions. This then allowed the analysis erf the data in die full three dimensional configuration where 
comparisons could be made in the multi-axis performance of die hardware. This chapter will 
describe the Leaijet flight experimental environment and die acceleration performance erf the isolated 
payload with respect to its dynamic inputs while attached to the aircraft during die low gravity 
maneuver. 


4.2.1 Aircraft Experimental Environment 

The Leaijet demonstration hardware was housed in a standard Leaijet rack. These racks 
have standard instrumentation interfaces with T-rail mountings to attach to the aircraft fuselage. The 
Leaijet hardware consisted erf two instrumentation racks, one for the levitated test section and the 
second housing die control computer, do power supplies, and support electronics. The levitated test 
section was interfaced with a trunnion support package housed internal to a standard rack allowing 



105 


the experimental package to pivot about a trunnion support shaft Figure 4.8 is a photograph of the 
two Learjet mounted experimental racks where die tnmnioned support is shown pivoted about its 
support shaft The electronics rack shows the control computer with the proximitor, accelerometer, 
and magnet dc power supplies. A master data acquisition computer was also housed in this rack 
which poled fourteen Learjet referenced data channels. In addition, a monitor and two current meters 
were attached to the electronics rack, where the current meters gave the total magnetic actuators ± 
current draw. The trunnion attached hardware consisted of die levitated platform, three actuation 
pods, the control sensors, and the magnetic actuators current control power amplifiers. Figure 4.9 is 
a photograph showing an end view of die tnmnioned payload. The top part of die tnmnioned cube 
housed the twelve power amplifiers and the proximitor drive signal conditioning circuits. The bottom 
of the tnmnioned volume housed the actual isolation system. The isolated payload consisted of a 
ferromagnetic structure where an autonomous six channel data acquisition system (DAS) was housed 
and slaved to the master DAS. The autonomous slaved system was time synchronized with the 
master DAS to take the acceleration time histories by triggering the two data acquisition systems 
using a physical trigger which was pressed prior to entering the low gravity portion of die Keplarien 
trajectories. 

In order to describe die experimental acceleration environment for the Learjet tests, the 
dynamic rigid body control of the aircraft will be discussed. The Learjet low gravity guidance system, 
provides low frequency acceleration data for the aircraft's X, Y, and Z axes as shown in Figure 4.10 
[23]. Acceleration levels are displayed in die cockpit for crew guidance and are recorded far research 
data correlation. However, attached experiments are subjected to higher frequency vibratory 
accelerations emanating from structure-home noise and direct acoustic energy. Additionally, since 
the reference accelerometers used for guidance are usually located near the aircraft's center of gravity 
(C.G.), experiments located at other locations experience different rigid body accelerations, especially 
in the Z axis, due to the rotation of the fuselage about the Y axis. 
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Figure 4.8: Learjet demonstration hardware. 





Figure 4.9: Leaijet trunniooed rack support structure and isolated platform. 
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Figure 4.10: Aircraft coordinate system. 



Figure 4.11: Aircraft rigid body dynamics coordinate systems definition. 




109 


The system by which these trajectories are flown may best be described by theoretically 
deriving the equation of motion for die flight maneuver in two dimensions, or assuming simple 
planar motion. A typical trajectory is illustrated in Figure 4.11 where the low gravity portion of this 
trajectory begins after a high g maneuver, as indicated in the figure. The coordinates far the 
trajectory and body fixed axes will be defined as diown in figure 4.11. The vector r is defined as the 
vector tracking the position of the trajectory reference sensor to an inertial coordinate system. In die 
body fixed frame, the distance from an accelerometer to this coordinate system is defined by a vector 
length t in the ej direction. Summarizing, the two acceleration equations for the body fixed 
coordinate system can be written as [24]: 

r, = (v£ - g„ cos qje r + ( v - sin Qe, (4.1) 

and, 

r, = (4 -t\-g B cos l)e r +(v -*4 2 ~ & sin ^)e , . (4.2) 

The rotational velocity or pitch rate and forward thrust of the aircraft is being controlled in order to 
cancel the acceleration due to gravity in the body fixed axes. This flight control is accomplished by 
viewing a Light Emitting Diode (LED) display where one axis of the display is proportional to 
acceleration in the radial direction and the other axis is proportional to tangential acceleration. Both 
of these accelerations are measured at the remotely mounted reference sensor. The maneuver 'control 
loop" is closed by the pilots controlling the pitch rate and thrust of the aircraft. 

Controlling the aircraft trajectory in this manner, the reference sensor, in effect, becomes the 
electronic center of gravity of the aircraft. As seen in equation (4.2), there is an additional 
acceleration term due to the pitch rate and the pitch rate as a function of time in the radial and 
tangential coordinate directions, respectively. These acceleration components are linearly related to 
the distance from a point to the electronic C.G.. This simplified planar equation of morion shows that 
the residual acceleration of an experiment can be substantially different from that of the control 


sensor accelerations. 
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In order to adequately describe the dynamics of a low gravity maneuver, and thus die 
experimental environment for the Learjet hardware, it is advantageous to give the filtered time trace 
of the accelerations throughout the aircraft as well as the acceleration spectr um during the off loaded 
phase of the trajectory. The filtered acceleration traces give die residual trends of a typical low 
gravity flight maneuver as well as the sensors acceleration history based on its attachment to the 
aircraft. In addition, these filtered time histories are used k> defined the three phases of a low gravfty 
trajectory. Figure 4.12, shows die filtered time history traces of three vertically positioned 
accelerometers whose relative distance, 2Z to 4Z and 3Z to 4Z, were a pprox im ately 0.838 m (33 in.) 
and 0.94 m (37 in.), respectively. The average mean acceleration in non-dimensional units for 
accelerometers 2Z, 3Z, and 4Z, were -0.0025, -0.0027, and -0,0021 g’s, respectively. The residual 
relative offsets shown in Figure 4.12 can be attributed to two sources, one being rigid body angular 
accelerations of the aircraft and the second due to non-calibrated Idas voltages inherent to the signal 
conditioning circuits. A correction for die later was attempted by zeroing the signal conditioning 
boards and taking an offset reading in addition to die bias calibration factors used for die QA-2000 
proof mass accelerometers. The increased negative residual acceleration offset with increasing 
distance from the control sensor location for sensors 2Z, 3Z, and 4Z, is consistent with die expected 
magnitude increase in residual average accelerations because of their location relative to die control 
sensor. To attempt a quantitative check of the measured observation, the second derivative of a pitch 
gyroscope time history was taken for die same time history giving a peak value of 0.3 (deg/s^). This 
ralnilatfd angular acceleration compared well with observed measurements of 0.27 and 036 (deg/s^). 
These measured values are calculated from die difference in mean acceleration values of sensors 2Z to 
4Z, and 3Z to 4Z, respectively. The calculations using die gyroscopic data assume only planar 
motion. Out-of-plane dynamics due Id roll and yaw were not accounted for. These residual 
accelerations were on the order of 10 to 20 percent of die residual g experienced at die control sensor 
for a typical accelerometer position of 234 m (100 in.) from the aircraft trajectory control point 
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A typical acceleration envir onment for the Leaijet tests will be defined by tile filtered time 
trace as shown in Figure 4.12. As shown on tins figure, the aircraft accelerations are shown having 
transitioned from the high g mane uver into the entry phase of die trajectory. This entry phase is 
defined by a posh over maneuver by the pilots, to zero out the pitch axis accelerometer. The aircraft 
typically overshoots the target "zero g" and an oscillatory residual acceleration condition follows until 
the pilot controlling the pitch axis of die aircraft stabilizes the control point. Subsequently, a 10 to 15 
second stable, low gravity, residual acceleration environment is established. At die end of die stable 
low gravity time the aircraft is pitched approxim ately 45 to 50 degrees nose down, at which time die 
pilots pull out of the trajectory and prepare for die next maneuver. The low gravity trajectories can be 
dissected into three major phases, entry, low gravity, and exit. The low gravity trajectory phase, 
which defines the testing environment, typically lasts 10 to 15 seconds. This stationary part of die 
trajectory therefore, limits the low frequency bandwidth of any isolation system as well as the testing 
bandwidth to approximately a 0.1 to 0.3 Hz low frequency boundary. 

The methodology used to produce the power spectral densities for each accelerometer during 
the low gravity trajectory first involved determining the times at which the stable low gravity period 
began and ended, determining bow many ensembles (time samples) could be processed into a 
spectrum average for that duration, and then processing die autospectra. A basic check of random 
process stationarity was conducted. It was demonstrated that there was minimal variation of the 
resultant power spectral densities if three or more low gravity ensembles were used. It was also noted 
that there was si gnificant variation in the power spectral densities if die entire low gravity maneuver 
acceleration data set were processed as a set of ensembles. This demonstrated the nonstationarity of 
the entry/low gravity/exit time domain process. An example of a longitudinal (X axis) 
accelerometer's time domain response, can be seen in Figure 4.13. Figure 4.13 clearly shows die drop 
in response level occurring during the low gravity period compared to die entry and exit maneuvers. 




Figure 4.13: X axis 
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The frequency spectrum magnitudes seen by the accelerometers of attached payloads are 
dependent upon the payload structural dynamics, as well as the aircraft's dynamic response during the 
low gravity maneuver. As demonstrated in figure 4.13, the aircraft's response during the low gravity 
phase of the trajectory is representative of an off-loaded structure, figure 4.14, shows the power 
spectrum of an accelerometer on the trunnioned experimental rack and one mounted on the actively 
isolated payload, calculated for the same flight as the residual traces shown in figure 4.12, for a 22 
and 12 second time history. The 22 second time history included both initial transient acceleration 
oscillations, and the low gravity portion of the trajectory, once the transient entry dynamics st a biliz ed. 
The 12 second time history only captured the quiet low gravity portion of the flight with the entry 
level accelerations removed. The time points per ensemble were set with a 50 percent Hanning 
window giving a minimum of three averages per power spectrum. This resulted in a frequency 
resolution of 0.14 and 0.17 Hz for the 22 and 12 second time histories, respectively. 

As seen in figure 4.14, the accelerometers attached to the experiment rack demonstrated 
responses both broadband in nature with indications of modal response specific to the trunnioned 
rack. As can be seen for the 22 second data set processed, the entry transient dynamics of the initial 
entry phase drtmmares the low frequency spectral response, while the 12 second data set shows a more 
broadband response with rack resonance's, for the rack mounted accelerometer. The isolation 
performance spectrums will be calculated and used for the transfer function calculations during the 
stable low gravity phase of the trajectories since the control system was not activated until the aircraft 
was into the offloaded environment. 

Inconsistencies in the fli ght regime and flight conditions such as aircraft weight caused 
differing structural responses, the altitude effects on engine output performance (RPM variations), 
and even the variability on the entry technique levels and styles from pilot to pilot will tend to affect 
the frequency vibro- acoustic response of the aircraft Even so, the individual channel responses 
demonstrated a degree of consistency from one trajectory to the next and allowed the testing of the 
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active isolation hardware in the full six degrees of freedom for the offloaded condition. Even though 
the low frequency bandwidth was limited due to the duration of die low gravity phase of the 
trajectories there was enough time to take statistically representative acceleration data in the sub- 
Hertz frequency regime. 

4 2.2 Aircraft Data Collection and Transfer Function Measurements 

The experimental data acquired during die active Learjet demonstration tests were obtained 
during a series of low gravity flight tests. These data sets represent a wide range of dynamic 
conditions which posed a unique challenge in acquiring low frequency acceleration data from a fully 
active six degree of freedom initially controlled payload. These flights were successful in 
demonstrating the sub-Hertz control of a payload in a uniquely dynamic oon-stationary acceleration 
environment The goal of isolating a payload in an orbital flight environment will be a much smaller 
task than the gathering of data during low gravity aircraft trajectories. 

The Learjet aircraft, trunnion frame, and die active isolation platform were instrumented 
with Sunstrand QA-2000 accelerometers, hi addition, a Humphrey V624-0825-1 two axis vertical 
displacement gyroscope was installed on the trunnion system to measure die rotation and pitch angles 
of the trunnion support structure with respect to the vertical, based on a ground coordinate system. 
The gyroscope had a resolution of ± 3 degrees (0.0524 radians). The data collection instrumentation 
is shown in Figure 4.15. Figure 4.15 illustrates the input acceleration sensor locations for die transfer 
function nalmlarinnc where the output data channels acceleration sensors were mounted on the 
isolation platform in the configuration shown in Figure 4.16. Each accelerometer block in Figure 
4.15 consisted of three accelerometers oriented othogonally and aligned with the three aircraft-fixed 
coordinates X, Y, and Z. As previously stated, a total of 20 channels of data were acquired during 
each trajectory. The data was stored with 14 bit digital resolution and a sampling frequency of 142 
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Figure 4.16: Output power spectrum accelerometer locations on platform. 
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Hz per channel. The acceleration measurement system had a 1 5 Volt dynamic range with a typical 
sensor sensitivity of S Volts/.5 g. Each accelerometers sensitivity was used in calmhumg the per 
channel sensitivity values in Volts per g. The resolution of this data acquisition sy s t e m was 61 pg per 
bit using a typical accelerometer sensitivity of 5 Volts per g. As previously mentioned, the data 
acquis i ti o n system was manually triggered during die on-set of die low gravity portion of each 
trajectory. The manual trigger was pressed when the operator felt the transition from die 25 g pull- 
up into the low gravity portion of die trajectory, at which time the data acquisition system was 
initiated and data was taken for 27.6 seconds. Manual triggering, die non-stationary trajectory, and 
trajectory acceleration histories, all make comparisons between trajectories diffic ult Therefore, 
comparisons must be made with reference to specific phases of a maneuver such as the stahilirwi low 
gravity portion in which die active system was activated. 

The low gravity parabolic maneuver is depicted in Figure 4.17, where each phase of the 
trajectory is shown, entry/ low gravity/ exit Referring to this figure, die active demonstration 
hardware was activated during the low gravity phase of the trajectory, stabilized, and data was 
acquired to calculate die frequency response of the payload, hi order to best recreate the dynamics of 
both the actively controlled payload and its support structure, two data acquisition systems were 
flown. A slaved autonomous six channel data acquisition system (DAS) was attached to die 
suspended platform, while a master 14 channel DAS was flown for die Leaijet acceleration and 
rotational environment histories. A total of 18 acceleration and two gyroscopic data chann els where 
digitized by the 14 bit converter. 

The natural frequency of the demonstration hardware was set at about 0.6 Hz for the 
following pitch axis transfer function comparisons. However, die demonstration flights posed 
considerable environmental challenges. Therefore, the system was intentionally over damped in 
order to insure the stabilization of the platform after die initial conditions during die push over {base 
of the low gravity trajectories. Figure 4.18 shows the frequency response curves for two typical 
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trajectories where the active system is under a closed loop, relative feedback and relative/inertial 
feedback/feedforward control, respectively. These frequency response curves are given for the vertical 
direction where the acceleration spectrum of die payload is compared to that of die support structure. 

The relative and inertial frequency response curves were calculated from 17 and 14 second 
low gravity time histones. In order to get fairly representative frequency response functions for both 
cases, the elements per ensemble, with a SO percent Harming window, were set to generate die {dotted 
curves with stable results. This gave a frequency resolution of 0.2 and 0.24 Hz for the relative and 
inertial cases, respectively. Figure 4.18 (a) through (d), show die response functions of a soft well 
damped system with a natural frequency of about 0.6 Hz. The inertially referenced curves as 
compared to the relative feedback curves show the system’s increased roll off and a tt e nuati on as a 
function of frequency. The expected increase in attenuation of inertial feedforward compared to 
relative control was masked in the bandwidth from 2 to 10 Hz for these sensor to sensor comparisons. 

The multi-dimensional transfer function calculations for an inertially referenced case and a 
relative feedback case will be given. The cross-axis calculations are not presented in that the input 
power spectrums for two aircraft fixed locations were used, where each location consisted erf two 
triaxial acceleration measurements giving the three translational components in the aircraft fixed X, 
Y, and Z axes, respectively. Two rotational input power spectrums were estimated for die trunnion 
and aircraft fuselage fixed locations by differencing the two accelerometers sensitive to the pitch and 
yaw motions of the aircraft. The third rotational component was estimated by cross differencing the 
two locations, trunnion fixed and aircraft fixed, since the roll sensitive accelerometers for each 
individual location were aligned along the same accelerometer defined axis. Referring to Figure 4.15, 
the input power spectrums in the aircraft fixed translational axes were simply used directly as input 
power spectral reference channels for the X, Y, and Z axis frequency response calculati ons. This 
allowed the comparison between attenuation responses of the isolated payload with respect to each 
axis at the four locations, where two locations are paired as aircraft fixed and the other as trunnion 
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fixed accelerometers. This gave a good estimate of die overall attenuation from the aircraft 
attachment input acceleration power spectrum to the payload response as well as the tr unni on support 
structure to payload attenuation. Therefore, locations 1 through 4 X, Y, and Z accelerometers were 
used as a single axis reference to the X,Y, and Z axis output power spectrums of the payload. These 
payload output power spectral densities were e»lralat«t by taking the six inertial data sensors on the 
payload, and from their geometrical placements, resolving the appropriate six inertial degrees of 
freedom of the platform. Referring to Figure 4.16, die following equations rep resent the accelerations 
that the individual platform accelerometers sense along their sensitive axes. 


a* = — (-Xcos30+rsm30+M) 

So 

(4.3) 

So 

(4.4) 

a, 7 = — ( X cos 30 + 7 sin 30 - p x /?) 

So 

(4-5) 

0,8 = —(2+0.,/? cos 30-<{>ji?sin3o) 

So 

(4.6) 

(z-M) 

So 

(4.7) 

Ojq = — (z + 0 x /? cos 30 + j/? sin 30) 
So 

(4.8) 


Using these equations the following non-dimensional translational accelerations can be resolved for 
the platform inertial accelerations in the X, Y, and Z directions as well as the rotational accelerations 
of the platform, where these rotational accelerations are in units of radians/sec^. (Note: go is the 
acceleration of gravity and R is the radial distance from the center of die platform, R=5 inches.) 
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So 



(4-9) 
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(4.13) 


(4.14) 


Therefore, by using equations (4.9) through (4.14) die inertial six degrees of freedom were calculated 


from the respective acceleration channels aj 5 through a^Q. Referring to Figure 4.15, the rotational 
mpnt power spectrums were calculated by using difference relations between the aircraft and trunnion 
fixed accelerometers as follows: 



(4-15) 

®a£ren^ “ 2^ ^ °io)* 

(4.16) 

U-l=-fp( a 6-‘hl 

(4.17) 

^=■^( 03 - 01 °). 

(4.18) 

fittnuumm ~ ^ q 

(4.19) 


pjairmtf _ 35* 5" 


(4.20) 



124 


Once these acceleration time histones were calculated from the respective individual 
acceleration sensor channels the stable low gravity time was determined from the individual time 
histories. This low gravity time slice of die trajectory is the acceleration time history which must be 
processed in order to calculate the individual autospectra and the respective single axis frequency 
response curves for die inertial and relative control data sets. 

The following two data sets were taken from the same flight test The first data set processed 
was under relative control where the second data set was processed with the same relative parameters 
as the first but the inertial control loops closed during die low gravity trajectory. Again, since each 
trajectory can have substantially different initial conditions, the data was processed by determining 
the stable low gravity portion of each trajectory, and at which time die active feedback and or 
feedback/feedforward loops were closed, during these trajectories. This was accomplished by using 
the acceleration time histories of the isolated payload and determining when die initial acceleration 
spike occurred, giving the time reference of the controlled isolation period. In addition to the 
beginning of the processible acceleration thm» history, die duration of tins acceleration data set must 
be determined because each trajectory is subject to different operating environments which therefore, 
determine the overall time of each trajectory as well as the controllability of each maneuver by die 
pilots. The termination of each acceleration rime. history was determined by reviewing tbe total 
acceleration time trace and viewing a second acceleration spike on die isolation system which noted a 
wall-to-wall contact Each acceleration channel for die respective trajectory was reviewed and the 
rime slice with no acceleration spikes for all six acceleration output channels was used to process die 
frequency response curves and their respective antospectra. This of course ment that the overall 
spectral resolution would not be the same for the two trajectories compared. However, if three or 
more averages were available die frequency response functions were stable. Since the frequency 
domain is used far these calculations the trajectories can be compared to one another in the 
bandwidth defined by die total time history processed and the sampling frequency of die data 
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acquisition system. Figure 4.19 (a) and (b) gives die overall acc el e r ation t un e trace for a vertically 
mounted accelerometer on die trunnion support frame, where die entry/ low gravity/ exit phases can 
be seen for die relative control and die inertial control trajectories, respectively. Figure 4.20 (a) and 
(b) show the ycelerati nn time trace of a vertically mounted accelerometer attache d to die isolated 
system showing the period of control as well as the end of die process time history to be used for both 
trajectories, respectively. These two time slices for die relative and inertial control were used to 
generate the following frequency response curves for the X, Y, and Z dimensions and their respective 
rotations. For brevity, only the frequency response curves for the locations and the rotational 
definitions, as defined previously, will be given in this chapter. The associated autospectral densities 
of each frequency response curve is given in Appendix F. 

The time history processed for the relative case had a total time window of 16 seconds where 
760 points per ensemble average were used giving a spectral resolution of 381 lines for a frequency 
window from 0.01 to 75 Hz. The inertial time history had a total window time of 9.5 seconds where 
600 points per ensemble average was used giving a spectral resolution of 301 lines for the same 
frequency window. F^rh frequency response curve was calculated with a 50 percent Hanning window 
and the relative data set used five spectral averages where the inertial data set had three spectral 
averages. This caused die inertial frequency curves to be slightly more choppy however, comparisons 
in the frequ enc y response curves could be marie over die ope rational bandwidth of 0.1 to 75 Hz. 
From the t otal rime history of 16 and 9.5 seconds for the relative and inertial control trajectories, 
respectively, a samplin g rate per channel of 142 Hz gives a total of 2272 and 1349 points available for 
a spectral average. Therefore, for the relative case with a 760 point per ensemble average, die 
frequency resolution limit is 0.19 Hz while for the inertial trajectory a 600 point per ensemble average 
gives a frequency resolution limit of 0.24 Hz. Figures 4.21 through 4.24 (a) and (b), give the 
frequency response curves for die X direction for die single input reference locations one through four 
for both trajectories, respectively. Figures 4.25 through 4.28 (a) and (b) give die frequency response 
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(a) Payload under relative control. 

Figure 4.19: Acceleration time trace of full trajectory for aircraft mounted accelerometer. 
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(b) Payload under inertial control 

Figure 4.19: Acceleration time trace of fuD trajectory for aircraft mounted accelerometer. 
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(a) Payload under relative control 

Figure 4.20: Acceleration time trace of fall trajectory for platform mounted accelerometer 
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(b) Payload under inertial control. 

Figure 4.20: Acceleration time trace of full trajectory for platform mounted accelerometer. 
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curves for the Y direction, and Figures 4.29 through 432 (a) and (b) give die frequency reh ouse 
curves for the Z direction, for the single input reference locations one through four for both 
trajectories. Figures 433 and 434 (a) and (b) give the frequency response curves for the 0* rotation 
for the trunnion location and the aircraft location, respectively. Figures 435 and 436 (a) and (b) give 
die frequency response curves far the rotations as defined in equations (4.17) and (4.18), and 
Figures 4.37 and 4.38 (a) and (b) are the frequency response curves for the ^ rotation as defined in 
equations (4.19) and (430). The response functions have been plotted with the same minimum and 
maximum values for both trajectories. This data will be discussed in die next chapter imder data 
analysis and discussions. 
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(a) Relative control of platform. 

Figure 421 : Transfer function calculation for platform x vs. aircraft location 1 x. 
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(b) Inertial control of platform. 

Figure 4.21: Transfer function calculation for platform z vs. aircraft location 1 x. 
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(a) Relative control of platform. 

Figure 422 : Transfer function calculation for platform x vs. aircraft location 2 x. 
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(b) Inertial control of platform. 

Figure 4.22: Transfer function calculation for platform x vs. aircraft location 2 x. 
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(a) Relative control of platform. 

Figure 4.23: Transfer function calculation for platform x vs. aircraft location 3 x. 
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(b) Inertial control of platform. 

Figure 4.23: Transfer function calcHlation for platform x to. aircraft location 3 x. 
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(a) Relative control of platform. 

Figure 4.24: Transfer function calculation for platform x vs. aircraft location 4 x. 
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(b) Inertial control of platform. 

Figure 4.24: Transfer function calculation for platform x vs. aircraft location 4 x. 
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(a) Relative control of platform. 

Figure 4.25: Transfer function calculation for platform y vs. aircraft location 1 y. 
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(b) Inertial control of platform. 

Figure 4.25: Transfer function calculation for platform y vs. aircraft location 1 y. 
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(a) Relative control of platform. 

Figure 4.26: Transfer function calculation for platform y vs. aircraft location 2 j. 
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(b) Inertial control of platform. 

Figure 4J6: Transfer function calculation for platform y vs. aircraft location 2 y. 
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(a) Relative control of platform. 

Figure 4.27: Transfer function calculation for platform y vs. aircraft location 3 y. 
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(b) Inertial control of platform. 

Figure 4.27: Transfer function calculation for platform y vs. aircraft location 3 y. 
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(a) Relative control of platform. 

Figure 4.28: Transfer function calculation for platform y vs. aircraft location 4 y. 
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(b) Inertial control of platform. 

Figure 4.28: Transfer function calculation for platform y vs. aircraft location 4 y. 
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(a) Relative control of platform. 

Figure 4.29: Transfer function calcnlation for platform z vs. aircraft location 1 z. 
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(b) Inertial control of platform. 

Figure 4.29: Transfer function calculation for platform z vs. aircraft location 1 z. 
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(a) Relative control of platform. 

Figure 4 JO: Transfer function calculation for platform z vs. aircraft location 2 z. 
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(b) Inertial control of piatform. 

Figure 4 JO: Transfer function calculation for platform z vs. aircraft location 2 z. 
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(a) Relative control of platform. 

Figure 4.31: Transfer function calculati on for platform z vs. aircraft location 3 z. 
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(b) Inertial control of platform. 

Figure 4J1: Transfer function calculation for platform z vs. aircraft location 3 z. 
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(a) Relative control of platform. 

Figure 4.32: Transfer function calculation for platform z vs. aircraft location 4 z. 
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(b) Inertial control of platform. 

Figure 432 : Transfer function calculation for platform z vs. aircraft location 4 z. 
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(a) Relative control of platform. 

Figure 4.33: Transfer function calculation for platform 0 Z vs. trunnion 6 r 
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(b) Inertial control of platform. 

Figure 433: Transfer function calculation for platform 0 X ys. trunnion 0 r 
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(a) Relative control of platform. 

Figure 434: Transfer function calculation for platform 6 X vs. aircraft 0 r 
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(b) Inertial control of platform. 

Figure 434: Transfer function calculation for platform 8 X vs. aircraft 0 T 
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(a) Relative control of platform. 

Figure 4.35: Transfer function calculation for platform ty vs. ty2*3* 
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(b) Inertial control of platform. 

Figure 4.35: Transfer function calculation for platform ty ^ ♦y2-3- 
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(a) Relative control of platform. 

Figure 436: Transfer function calculation for platform ty +yl-4« 
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(b) Inertial control of platform. 

Figure 4.36: Transfer function calculation for platform ^ vs. tyl-4< 
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(a) Relative control of platform. 

Figure 437: Transfer function calculation for platform p z vs. trunnion p z . 
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(b) Inertial control of platform. 

Figure 4.37: Transfer function calculation for platform P z vs. trunnion p r 
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(a) Relative control of platform. 

Figure 438: Transfer function calculation for platform p z vs. aircraft P 2 . 
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(b) Inertial control of platform. 

Figure 438: Transfer function calculation for platform P z vs. aircraft |3 Z . 
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V. Development and Demonstration Systems 
Experimental Measurements Analysis 


Two six degree of freedom active magnetic isolation systems were designed and developed 
for die demonstration of inertial sub-Hertz isolation. The laboratory hardware was de sign ed to have a 
configurable system where multiple control algorithms and inertial isolation techniques could be 
tested in a ground based environment The proto-type hardware was tested in the h o ri zo n tal plane 
due to constraints associated with tire ground based 1 g environment The Learjet demonstration 
hardware was designed as a functional equivalent to die laboratory system demonstrating multi-axis 
control and testing in an off loaded environment during low gravity Keplarien trajectories. The 
following sections will cover the experimental performance measurements data analysis, and 
discussions of frequency response trends, and system acceleration and attenuation performance limits. 

5.1 Laboratory Hardware Acceleration Measurements and 
Transfer Function Analyses 

The noise floor acceleration measurements for die payload fixed laboratory accelerometers 
were shown in Figures 32 and 3.3. These curves show die horizontal x and y accelerations of the 
payload where the space station "micorgravity" requirement was plotted for reference. Estimates of 
system noise for the laboratory hardware predictions of the total system noise performance in g's 
verses frequency can be made. Referring to Chapter 3 section 3.1.3, die noise floor for relative 
position control was given, where die eddy-current probes signal conditioning circuits had a noise 
floor performance of approximately 10"^ to 10~6 volts nns over the bandwidth from 0.1 to 100 Hz. 
Since the frequency response of the eddy-current probes over the control bandwidth is flat, with the 
signal conditioning filters having a gain of approximately 4.0, the sensitivity is estimated at 80 V/in 
from a nominal transducer sensitivity of 20 V/in. The horizontal 12 bit digital converters with a ± 5 
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Volt dynamic range gave a voltage per count sensitivity of 10/2*2 v/coanL Therefore, die eddy- 
current and control system deadband can be estimated at 3.05x10^ inches. This is equivalent to a 
voltage noise of the controller of 2.4x10*3 Volts peak-to-peak. In comparison to the eddy-current 
probes voltage, noise becomes die limiting factor for die steady state case. Using this dimens i on al 
resolution of 3.05x10*3 inches, one can estimate an acceleration performance for steady state 
operation. Taking die riftarfhmri of the controller both steady state and in die controller bandwidth, a 
deadband resolution of 3.05x10*3 inches is an estimate of die amplitude peak to peak resolution 
across the controller bandwidth. If the control loop was to command die dimensional position of die 
payload across its bandwidth to die deadband resolution of die digital system, a non- dimens ional 
acceleration verses frequency curve can be calculated using die deadband resolution esti ma te, where 
a nmit (B>) = (3.05x10*3 inchesK27rf bandwidlh ) 2 /g 0 . Figure 5.1 gives this curve with the acceleration 
measurements of the platforms x and y dimensions. The estimated noise limit only predicts die sub- 
Hertz response of die acceleration noise floor measurements. This acceleration limit es timate 
assumes there is a count change at die frequencies in the control bandwidth and leads to a 
conservative estimate for the deadband response, with no input excitation, as demonstrated by the 
support acceleration measurements in Figure 5.1. If one assumes the signal conditioning circuits and 
the analog-to-digital converters are stable, then the digital-to-analog converters output will dictate the 
dimensional stability of the active system with no feedback signal. Using (he displacement probes 
sensitivity of 80 V/in and die voltage noise from the digital-to-analog converters, a dimensional 
stability can be estimated as shown in Figure 3.18. Using this spectrum and taking the second 
derivative in the frequency domain, one arrives at die acceleration noise spectrum shown in Figure 
52 . This noise measurement was done on die 16 bit hardware, therefore, if an adjustment of 2* were 
marl** on this curve there is good agreement with the acceleration measurements shown in Figure 5.1 
for the frequency range from 2 to 10 Hz. The noise floor estimate using the input dimensional 


5/3/91 9:42AM 

SloH=O.OMi Spon=l?.500Hi lint$=40l No (Mp Honninq 


169 



Flgore 5.1: 


Acceleration noise performance spectra for laboratory system with 
performance boundary estimates. 
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resolution is an es timate of the limiting case at the control systems set natural frequency or the 
controller steady state resolution. The systems proportional feedback term controls the system at its 
set natural frequency to the dis pl acement resolution of the feedback sensors. The peak-to-peak 
response of die 12 bit system with a natural frequency of 0.65 Hz can be estimated as follows: 

Om. = (27K0,) 2 . (5.1) 

So 

Using the 12 bit dimensional deadband and die systems natural frequency an estimated one count 
change in the controller gives an acceleration peak-to-peak stability of 1317x10-6 g’s. However, the 
controller must see a few count changes for control of the platform giving a value of 2-633x10-6 g’s 
where this estimated acceleration stability is in good agreement with die amplitude response of the 
measured values in Figure 5.1 at the systems natural frequency. 

The transfer function results for the prototype hardware are dictated by the digital resolution 
of the input converters. Referring to the transfer function calculations from the input power and 
output power spectral densities for the horizontal plane, shown in Figure 4.7, the attenuation profiles 
for the inertial control of the platform demonstrated a 100 dB/decade roll-off with a stop frequency of 
about 2 Hz. The attenuation bottomed out at approximately -30 dB where the transfer functions 
leveled out This leveling out of the transfer function can be explained using die dimmsinnal 
resolution of the control system by analyzing die input and output power spectrums used in the 
transfer function calculations of Figure 4.7. The input and output spectral densities from the martial 
transfer function calculations are {dotted in Figure 53. This figure shows the output spectrums being 
attenuated where they bottom out and follow the dimensional resolution of the controller plotted on 
the same figure. The acceleration feedforward resolution for the laboratory hardware was 
approximately 6363xl0~^ g’s. Figure 53 shows the input spectrum crossing this feed fo rwar d 
resolution at approximately 05 to 1 Hz where the inertial cancellation begins to take effect and die 
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Figure S3: Input and output power spectra for inertial transfer function calculations. 
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output power spectrams begin to attenuate the input accelerations of the support structure within the 
dimensional resolution of the feedback sensors. This control approach, as m e nti o ned in the 
theoretical developments, commands the payload to be inertially still where an observer on the 
support structure would find the payload in motion equal to the resolution of the feedback oontroL 
Therefore, the steady state and controller bandwidth noise performance is bounded by die digital 
resolution of the control loops and by the output resolution of die digital-to-analog converters, 
assuming no input feedback signals. The transfer function is also bounded by the digital resolution 
where, theoretically, for an infinite resolution system and no sensor to sensor phase shifts, die control 
system should command an infinite attenuation or to die sensor noise floors in the 
feedback/feedforward control system. 

Figure 5.4 depicts the relative case with die input and output power spectrams plotted. The 
dimensional resolution is superimposed showing the output power spectrams following this curve as 
in Figure 5.3. The space station requirement is also plotted in Figures 5.3 and 5.4 as a reference. 
These sets of data show the advantages of inertial and active isolation in general for low frequency 
requirements. Specific to "microgravity" levels, it is quite apparent that although laboratory digital 
equipment was successfully utilized for the developments of these state of die art isolation techniques, 
at least 16 bit high resolution control is dictated by the stringent "microgravity" requirements for 
orbital platforms. 

5*2 Learjet Hardware Acceleration Measurements and 
Transfer Function Analyses 

The Leaijet demonstration hardware was designed as a functional equivalent to the 
laboratory prototype system. Both systems were actively controlled in six degrees of freedom where 
the demonstration hardware was built to investigate the off loaded response by testing during 
Keplarien low gravity trajectories. This environment posed unique limitations in testing however, it 
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was adequate for demonstrating the sub-Hertz inertial control of a package, where the final v alidation 
of such technology to die 0.01 Hz range will be accomplished during orbital flight 

The experimental data taken during the low gravity trajectories encompassed the full 
trajectory acceleration phases, entry/low gravity/exit, where the actively controlled payload was 
activated during the oscillatory part of the entry phase as depicted in the full acceleration trace of 
Figures 4.19 (a) and (b). The total current draw for the twelve magnets during th e higb-g portion of 
the trajectories was calculated at approximately 14 amps for a 23 g pull-up maneuver. Howe v er, the 
g's pulled during each trajectory vary from 2 to 3 g’s and therefore, the current draw to maintain die 
equilibrium gap in many instances was substantially greater than 14 amps which tripped the 
experimental power circuit breakers, where the max current draw allowed was 15 amps. Therefore, 
the control computer was constantly active with all input sensors powered. However, the power 
amplifiers were not activated until the pull-up portion of the entry phase was complete and the 
trajectory was into the free-fall condition thus, offloading the reference frame of the aircraft This of 
course gave an initial shock to the payload. Therefore, die control parameters were set to apply 
substantial damping to the system to maintain control and stabilize the platform in one to two seconds 
after power was applied. This then typically allowed a window of 10 to 15 seconds to acquire good 
actively controlled payload acceleration (fata far die calculation of the transfer functions plotted in 
Chapter 4. 

In order to acquire a quick estimate of the attenuation achieved, there were four transfer 
functions calculated for an inertial control case and a relative case. These four transfer functions 
were calculated by referring an isolated payload pitch axis accelerometer to the four aircraft fixed 
pitch axis locations. The four tocations were defined in Figures 4.15 (a) and (b). The pitch axis or Z 
axis accelerometers were used as a quick performance check. This axis was chosen doe to the largely 
non-stationary acceleration environment caused by the dynamics of the trajectory maneuver. The 
curves plotted do not demonstrate the rapid roll off achieved during the laboratory 
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measurements because of the large relative damping required to provide a stable entry phase and 
therefore, giving a long enough window of time to acquire data for calculating die sub-Hertz power 
spectral densities and their associated transfer functions. However, in demonstrating active inertial 
verses relative control the Learjet aircraft provided an ad e q nat e means for a qu antita tive evaluation of 
these active twhniqnfts in an off loaded case. Figures 18 (a) and (b) depict (be a t te nu ati o n of die 
platform pitch axis accelerations with respect to the direct support structure provided by die 
tnmnioned frame. These two plots show the platform attenuating the support structure input power 
by approximately -20 to -30 dB at 10 Hz. The inertial control curves demonstrate the inertial 
cancellation approach where at the lower frequencies the payload does not want to follow its support 
Therefore, the ntiennatian of the active system is dropping off well before die relatively referenced 
payload. The frequency range from 3 id 10 Hz is masked by direct vibro-acoustic disturbances where 
the payload will react similarly in die feedforward or feedback control cases. Figures 4.18 (c) and (d) 
are the transfer functions ral mlateri for die pitch axis accelerometer referenced to locations three and 
four directly mounted to die aircraft fuselage. The same response for die inertial and non-incrtial 
control cases are seen as in Figures 4.18 (a) and (b). 

The purpose of the Learjet demonstrations was to quantify the a ttenuat i o n performance of the 
active magnetic system in multi-dimensions for an off loaded case. This was performed by resolving 
die six inertial degrees of freedom of the isolated platform and calculating the transfer functions for 
these six spectrums referenced to the trunnion and aircraft fixed accelerometers. The translational 
a ccel erati ons for the aircraft and trunnion fixed locations were used for the XYZ transfer functions 
while die rotational accelerations were estimated by the geometrical differencing of die appropriate 
mpnt power acceleration sensors as defined in section 422 . Figures 4.19 and 4.20 (a) and (b) show 
the full time traces of the relative and inertial acceleration traces for locations one and six, as defined 
in Figures 4.1S and 4.16, which were used to compute the multi-dimensional transfer functions. 
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Performing the inertial control during the parabolic flight posed considerable difficulty since 
die total trajectory could not be controlled due to die large current draw needed during the high-g 
man euver. The inertial circuits were set with a bandpass beginning at 0.03 Hz, where for die Lear jet 
environment was a liability. This liability was found once die control actuators were powered and the 
inwtiai feedforward sensors sensed this large acceleration. The long time constant assoc i a t ed with 
die 0.03 Hz high pass filter required a few seconds to reach die steady state co n dition, limiting die 
isolation window. Therefore, die inertial loops were not closed for a few seconds in order to allow 
this t ransient condition to subside. Therefore, die controller would not co mmand die isolation system 
to react to this large input which was self induced by die Learjet oper ationa l procedure. Although 
this caused problems in the ultimate performance for die inertial control, comparisons could be made 
in the low frequency decoupling of the isolated payload with its support environ men t while the 
inertial loops were closed. The relative time histories processed had a 21% greater frequency 
resolution and allowed die transfer functions to be calculated with five averages while the inertial 
data permitted only three averages. 

Comparing the X axis response from Figures 4.21 through 4.24 (a) and (b), where the 
relative control cases are shown in Figures 4.21 dvough 4.24 (a), die inertial cancellation may be 
seen to command the platform to be decoupled to its support at low frequencies. However, the lack of 
ave rag es and the direct disturbances caused by airborne noise and the initial conditions mask the 
performance in the 2 to 10 Hz range because a data acquisition card cage exhibited a strong 
resonance at approximately nine Hz. In this case die relative loops are more effective at d amping out 
this direct excitation source. This response, however, is specific to the testing environment of the 
aircraft. The roll off of the inertial control is quantifiably greater as shown in the plots. The relative 
loops show a well damped low frequency system with an approximate -10 dB/decade slope past 3 Hz 
while the initial control with die same relative parameters exhibits a 30 to 40 dB/decade roll off past 
10 Hz with the low frequency response being a tt e nuat e d at -10 to -20 dB. 
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The Y axis response of the relative and inertial transfer functions from Figures 425 to 428 
shows the relative curves demonstrating a damped 2 Hz frequency filtering into the payload, whereas 
the inertial loops are canceling this response with a fairly flat -20 dB attenuation from 0 5 to 75 Hz. 
This 2 Hz spectrum is believed to be caused by die pilots reaction time of a ppro xim ately 0.5 to 1 
second while controlling die pitch axis of die airplane. 

The Z axis transfer function calculations for both cases {dotted in Figures 429 and 432 (a) 
and (b) for the relative and inertial control , respectively, exhibit die same overall response as the X 
and Y axes, where die low frequency attenuation is increased by die inertial loops. The br oadband 
attenuation of the inertial response shows the direct energy forcing die acceleration response of die 
payload where the relative control does a better job of damping die energy input to the payload. This 
of course was expected since only inertial feedforward control is effective for base disturbance 
cancellation. 

The rotational transfer functions for die 0 X , 4y> and degrees of freedom showed 
approximately the same response for both control cases. The differencing of the accelerometers to 
estimate the rotational accelerations is not highly accurate. However, for the P z and transfer 
functions the inertial control {dots did demonstrate the low frequency decoupling of the payload 
whereas the 2 to 10 Hz frequency band was madrwl by the direct disturbances. 

It was very difficult to predict and estimate die dynamic environment of die aircraft due to 
die non-stationary aspects of the maneuver as well as the non controlled inputs to this environment as 
was mentioned previously. However, quantitative comparisons were made assuming a fairly 
representative low gravity environment for both trajectories used in each comparison during the stable 
low gravity phase. Of course, initial conditions and exit conditions effect the power spectrums and 
with a limited window of active data a limited number of averages for die spectral densities mates the 
low frequency regime difficult to compare trajectories. 


VI. Conclusions 


In order to utilize the low gravity environment of space a well characterized and 
reproducible acceleration environment will dictate the isolation of many true "microgravity’ sensitive 
payloads. It has been found that many fundamental materials and fluids experiments are limited in 
their potential success due to excessive dynamic excitations present on die Shuttle orbiter and on any 
future or existing orbital carriers. Though the frequency content will be vehicle specific, all orbital 
vehicles are fairly soft structures and will be excited by manned presence as well as various dynamic 
pieces of support equipment such as pumps, motors and orbit keeping thrusters. The most difficult 
portion of the dynamic disturbances present on die various orbital vehicles are the sub-Hertz to 10 Hz 
range where, for the Shuttle, the crew accounts for many of these excitations while the future larger 
space structures may have many structural resonances below 1 He as welL 

The need for active verses passive isolation has been discussed thoroughly for the past 
several years far the orbital environment. Although, for certain specific cases a passive system may 
be sufficient for a certain bandwidth of interest, this work has been focused on die sub-Hertz or low 
frequency isolation problem where an active system is dictated. This work has focused on two phases 
of the "microgravity" vibration isolation research to date, development of theoretical approaches and 
demonstration of a fully active six degree of freedom system which illustrates the advantages of active 
inertial and relative control of a package. The inertial referencing of a payload has many advantages 
over relative feedback soft suspension systems. Therefore, an active six degree of freedom prototype 
development system was designed and built to verify theoretical predictions, as well as to provide an 
experimental testbed for die study of different controllers. This development system was then 
transitioned into a higher fidelity demonstration package for the testing of a feedback/feedforward 
fully suspended system during a reduced gravity environment The successful demonstration of such 
devices in a fully off loaded case was a good test of robustness, given the severe limitations posed by 
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the testing environment during low gravity parabolic flights. The full ve ri fi c atio n and technology 
dem onstration of these devices will have to occur during orbital flights, since this is die int e nd ed 
environment for these isolation platforms. Testing in die low frequency regime of 0.01 to 0.1 Hz is 
not practical in a ground based p rog ram which is another limitation dial can only be overcome by 
space based testing . The p ro t ot y pe hardware successfully demonstrated the advantages of inertial 
feedforward cancellation where a 100 dB/decade roll off was achieved. The relative system was also 
successfully demonstrated with it's inherent advantages over passive systems through die low 
stiffnesses achievable by decoupling of die static and dynamic stiffness of a system. The Learjet tests 
successfully demonstrated a billy active feedback and fe ed fo rward system in a low gravity 
environment and successfully reduced the excitation environment by at least one decade. The inertial 
loops closed during the Leaijet testing de mo nstrated die low frequency decoupling of the active 
system from its support environment However, in die Learjet vibro-acoustic environment an 
inertially damped system by inertial feedback would have been more effective in attenuating this 
specific environment 

6.2 Future Research 

Wi thin the inter national space community there is a considerable amount of research and 
development cm going relative to vibration isolation technology for microgravity payloads. NASA's 
Microgravity Sciences and Applications Division of OSSA sponsored a workshop to bring die various 
organizations involved in this area together. An International Wotkshop on Vibration Isolation 
Technology for Micxogravity Science Applications was held in April 1991, which included the 
mt(M-natinnai SSF partners, to review current technology status and assess future technology needs. 
All of the SSF international partners have active programs in this technology area. The European 
Space Agency (ESA) is currently developing hardware concepts through two contracts. The 
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hardware being developed is a Columbus facility for improving the microgravity environment quality 
for payloads. The facility is called the Microgravity Isolation Mount (MGIM) and uses current state 
of die art vibration isolation technology as a baseline. Validation in only one degree of freedom has 
been completed. ESA's work in this area has been published in die open literature. 

The Japanese Space Agency (NASDA) is also involved in vibration isolation technology 
development. The NASDA program is similar to existing NASA programs co ns i s ti n g of ground 
based testing and use of low gravity aircraft flights. NASDA however has not published much data 
on their efforts. The Canadian Space Agency (CSA) effort has been principally directed toward 
improving the low gravity environment aboard the KC-135 or any other aircraft with a Large Motion 
Isolation Mount (LMIM). There is merit in doing this since use of low gravity aircraft flights is a 
cost effective way to conduct noteworthy low gravity research. None of the on-going international 
programs are presently addressing the issue of flight evaluation and validation of vibration isolation 
technology and advanced concepts such as inertial referencing through feedforward and feedback 
means. At the International Workshop this issue was discussed with the conclusion that until this is 
done there will always be doubt about the suitability of orbital flight for science experiments and 
commercial processing sensitive to foe vibration environment An experiment has been submitted 
under an announcement of opportunity through foe Office of Advanced Concepts and Technology 
O ACT in the vibration isolation area and is focused on addressing this need. The ground based work 
that has been conducted in this area is extensive and foe results obtained indicate that foe state of foe 
art for vibration isolation can be extended into foe sub-Hertz regime. The ability to do so in a 
microgravity environment must be evaluated and verified in orbit. Successful veri fic ati o n and 
documentation of this capability will give users confidence to use orbital flight opportunities to 
attempt sensitive science endeavors and will give to commercial space initiatives the technology to 
provide a stable, well-characterized, affordable and reproducible microgravity environment critical to 
commercialize a reproducible product. The technology being evaluated will provide a data base 



enabling flight hardware developers to include in their hardware die specific components and control 
systems best suited for their requirements with confidence of success. 
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Appendix A 


Nomenclature 


a 

B 

*2 

c 

fc. 

F 

Fs 

Fo 

& 

G a 

G aifb 

sr 

G vifb 

G viff 

G lh-3h 

g 1v-3v 

g 

go 

h 

h lr-3r 

w 

i 

*b 

K 

f 

\ 


cross-sectional area of magnetic pole face 

average acceleration of gravity on die Earth's surface, 9.81 m/s^ 

peak amplitude 

servo force per payload mass 

inertial electromagnetic damping coefficient. Nsec/m 

relative electromagnetic (lamping coefficient. Nsec/m 

accelerometer voltage output, V 

direct disturbance, N 

isolator force, N 

force due to gravity on a system of mass m 

force exerted by magnetic pole face 

accelerometer's natural frequency 

acceleration feedback gain term in Appendix B 

acceleration feedback in text gain team 

inertial position feedforward gain term 

inertial velocity feedback gain tern in Appendix B 

inertial velocity feedback gain term in text 

inertial velocity feedforward gain term in text 

electro-magnetic actuator constant for horizontal actuators 

electro-magnetic actuator constant for vertical actuators 

non-dimensional acceleration 

acceleration due to Earth's gravitational field 

air gap between pole face and armature 

static equilibrium gap length 

relative horizontal displacement of horizontal actuation points one through three 

electromagnet current, a velocity 

current through a coil 

magnetic circuit current bias 

passive stiffness coefficient, NAn 

magnetic circuit current am plifier stiffness 

magnetic circuit isolator stiffness, NAn 

magnetic circuit position gain 

magnetic circuit current stiffness 

magnetic circuit sensor amplifier gain 

magnetic circuit velocity feed b ac k gain 

magnetic circuit position stiffness 


i length of pivoted beam 

lj actuator location from pivot point 

m mass, kg 

N number of ampere turns 

R resistance, SI 

T transmissibihty 
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f 


t 

u 


v lr-3r 

Xf(s) 


*r 

y r 

*r 

Pr 

e 

Or 

»o 

? 

§1 

“*1 


*2 

<D 


*F 

co 

“n 

<Ona 


time, sec 
position of base 

relative vertical displacements at vertical actuation points one through three 
mobility effectiveness 
position of payload 

relative displacement of payload in x direction 
relative displacement of payload in y direction 
relative displacement of payload in z direction 
relative angular displacement of payload about y axis 
angular displacement 

relative an g ular displacement of payload about x axis 

permeability of free space, 4nxl0*^ H/m 

control voltage 

proportional feedback voltage 

active damping coefficient 

passive damping coefficient 

time constant of sensing circuit 

time con stant of differentiator 

magnetic flux 

relative angular displacement of payload about z axis 

magnetic field strength 

excitation frequency 

system natural frequency 

active system resonance frequency 



Appendix B 


Theoretical Evaluation of Several Active Feedback and Feedforward Methods 

The following analysis su ppor ts the background summary of die inertial isolation techniques 
which, with the needed relative support information, lead to the prototype design. The approach 
taken for the response analysis on these one degree of freedom inertially based isolation systems were 
based cm weak by D. Schubert, Barry Controls Eastern Operation, Barry Wright Corporation, 
Watertown, Massachusetts. Three physical models. Figure B.l, will be analyzed to characterize 
isolation using acceleration, velocity, and both acceleration and velocity as the control feedback 
signals. 

For the physical systems described in Figure B.l two system disturbances exist The first is 
the base motion or a structural excitation, u. The second is an applied farce, F, to die payload mass, 
m. Both terms u and F are considered to be functions for which Laplace transformations can be 
realize d. The servo farce, F s , on the active control configurations is generated by sensing velocity, 
acceleration, or both and driving an actuator out of phase with the control signal. 

The sensor function is to convert die velocity, dx/dt (or acceleration, d^x/dt^), of the payload 
into an electrical voltage proportional to die excitation. The sensor output voltage, Ey, is proportional 
to velocity dx/dt (or acceleration, d^x/dt^) such that Eyy^Aydx/d t (or E va =A a d^x/dt). The voltage 
Ey is amplified with gain B v (or B^, such that the output voltage from the amplifier Egyy (or E^) 
will be Ej^^AyBydx/dt (or E ava =A a B a d^x/dt^), respectively. The voltage is then applied to a 
Lorentz force actuator coil of resistance R(Q) having negligible inductance. This coil is immersed in 
a magnetic field having a field strength, yf. Application of die voltage to the coil produces a current 
having magnitude I avv =E a yy/R (or I a y a =E ava /R). The effect of the current flowing through a coil of 
N turns of wire results in a force, F s , where Fgy= yNIgyy (or F^ V^^ava)- To simplify all of the 
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CD) 



(C) 


Figure BJ: 


(a) Velocity feedback. 

(b) Acceleration feedback. 

(c) Acceleration and velocity feedback. 

Physical description of three one degree of freedom inertial feedback 
configurations. 
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above, a gain term Gy (or G a ) is applied to the feedback signal. This gain term relates the velocity 
dx/dt (or acceleration, (dVdt 2 ) to die force (or F a ) such that F^sCydit/dt (or F^sG d^dt 2 ). 

Now FgyS^jydx/dt and F a = -G a d^x/dt^ , where die negative sign denotes negative feedback The 
sign applied to the servo force term, F r is governed by the direction die servo current flows thro ugh 
the coil of wire. 

The equations to model the three isolation systems are as follows. The active kniarinn 
system(s) differential equation of motion is obtained horn die force balance method wring Newton's 
first and second laws. Here the base motion, u, is actually a rime f unctio n, so o=u(t), with the gwne 
implied for the applied force, F, such that in actuality F^F(t). Thus, for velocity feedback. 


m^t = F+F„+K(u-x) 

(B.l) 

Substituting and rearranging give the following equations: 


d 2 x _ _ dx 
m-^ = F-G,- + m-x) 

(B2) 

m^-+G,—+Kx = F+Ku. 

dt 2 dt 

(B3) 

For acceleration feedback. 


m£± = F+F„+K(u-x) 

at 

(BA) 

Substituting and rearranging give the following equations: 


d 2 x _ _ d 2 x Trf 

m W =F 

(BS) 

d 2 x d 2 x 

m— r +G a ^+Kx = F+Ku 

dt 2 dt 2 

(B.6) 


For acceleration and velocity feedback, 

d 2 x 

m— T = F+F M2 +F„+K(u-x) 

dt 


(B.7) 
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Substituting and rearranging give the following equations: 

m ^ = F ~ G ‘^- G '% +KiU ~ X) M 

m^+G a 4-?-+G r —+fCx = F+Ku (B.9) 

dt 2 a dt 2 ' dt 

B.l Base-Excited Vibration Response 

To determine the base-excited response for die three vibration isolation systems, the force 
term F(t) is set to zero. The base displacement term u(t) is assumed to be a sine wave having a peak 

amp litude of B. Thus, u(t) = Bsin( t) where is a frequency tom. The three equations of motion 

become, for velocity feedback, 

m^-+G t —+Kx = KBsm((iX) (B.10) 

dt 2 dt 

for acceleration feedback, 

m^-+G a ^-+Kx = KBsm((Ot) (B.ll) 

dt 2 dt 2 

and for both acceleration and velocity feedback, 

m^-+G a ^-+G,—+Kx = KBsm((at) (B.12) 

dt 2 e dt 2 ' dt 

The equations are then transformed into the frequency domain using the Laplace 
transformation; 

m 

F{s) = \F{ty*dt (B.13) 

The velocity system equation becomes 

ms 2 X(s) + G r sX(s)+KX(s) = KU{s) 


(B.14) 
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the acceleration system equation becomes 

ms 2 X(s)+G a s 2 X(s)+KX(s) = KU(s) (B.15) 


and the acceleration and velocity system equation becomes 

ms 2 X(s) + G/X(s) + G^s) + £X(s) = KU(s) (B.16) 


The base-excited system transfer function is defined as X(s)/U(s) = T(s). Thus, die transfer 
functions for the three isolation systems are, for velocity feedback. 


Xjs ) _ K 

U(s) ms 2 +G r s+K 


(B.17) 


for acceleration feedback, 


X(j) _ K 
U(s) ms 2 +G a s 2 +K 


(B.18) 


for both acceleration and velocity feedback, 

K nl9) 

U(s) ms 2 + Gj 1 + Gj + K 

Transforming these transfer functions into a standard vibration notation gives, for velocity 
feedback. 


£/,)„ < 

U K j 2 +2 £cd ii $+cd 2 


for acceleration feedback. 



S 2 +-^ 5 - s * +<D - 


(B20) 


(B21) 


and for acceleration and velocity fee db a c k, 
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U 


(s) = - 


cat 


s 2 +-£rs 2 +2fyo M +(o 2 a 

SWo 


(B22) 


wberc Ga/m = Gyfa = 25 %, % = lf2CyOJKin) 1 ^. 

In this form the denominator of each t ransfe r function is called the characteristic equation. 
If the characteristic equation has real negative roots, then the vibration isolation system will not 
osc illate ; if it has complex roots, it will oscillate. If the base motion displacement time fraction u(t) 
is a "pure sinusoid", the steady state frequency response in complex form is given by letting s= jco. 
Thus, the frequency response is obtained by the function 

= 023) 

where s = jto and j = V~T. This frequency response, which is a vector, becomes, for velocity 
feedback. 


U J(S> ((D* - GO 2 ) + 2fyd m j(£> 


(B23) 


for acceleration feedback. 


u 


co: 


GO 


and for both velocity and acceleration feed ba c k. 




(B24) 


*(*>)= 

u 


co: 




+ 21,(0 Ja 


<B25) 
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The magnitudes of vibration measured on the isolated payload resulting from the shmatAfal 
excitation Bsin(ox) is die vector of length |X(yco)/£/(yO>)j. This value is a since the phase 

angle is not used and it is called the timsm usability. It is generally written as T(j<D) = 

|X(yo)) /£/(yco)j. Thus, the transmissibility function for die systems of interest become, for 

velocity feedback, 

1/2 

0*27) 

for acceleration feedback, 

1/2 

(B28) 

and for both acceleration and velocity feedback, 

1/2 

(B29) 

The data shown in Figures B.2-B.4 present transmissibility vs non-dimensional frequency as 
functions of the damping term, £ , and the acceleration term b 2 / 2 xfi 0 . 
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Figure B.2: Inertial velocity feedback transmissibility curves. 


196 



Figure BJ: Inertial acceleration feedback transmissibility curves. 
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Figure B.4: Inertial acceleration and velocity feedback traasm issibility corves. 
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B.2 Force-Excited Vibration Response 

Motion of the isolated payload can result from two excitation sources. The first is base 
motion. The second results from external forces applied directly to the isolated payload. Referring 
back to the first differential equations of motion (B3), (B.6), and (B.9), and setting the base 
excitation term u(t) to zero, the external force, F(t), is allowed to excite the payload. The equations of 
motion are, for velocity feedback, 

m^ r = F(t)-F„-Kx (B30) 

for acceleration feedback, 

m^ = F{t)-F m -Kx (B31) 

and for both velocity and acceleration feedback, 

m ~^ = F(t)-F K> -F„ -Kx (B32) 

These equations can be placed in the Laplace operation format, for velocity feedback, 

ms 2 X(s)+G v sX(s)+KX(s) = F(s) (B33) 

for acceleration feedback, 

ms 2 X(s ) + G^Xis) + KX(s) = F(s) (B34) 

and for both velocity and acceleration feedback, 

ms 2 X(s) + G a s 2 X{s) + G,sX(s) + KX(s) = F(s) (B35) 

The sensitivity of die isolated payload to the disturbing force, F(s), is characterized by a term 
called the isolated payload mobility. Mobility is die vector magnitude of X(s)/F(s). Mobility 
measures the amount by which die payload is deflected per unit of externally applied force. In 
T a place notation form, the equations for mobility become, for velocity feedback, 
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for acceleration feedback. 


X( v 1 

F KS,v ms 2 +G.S+K 


X ( . 1 

F (S) ‘ ms 2 +G t s 2 +K 


(B36) 


(B37) 


and for both acceleration and velocity feedback; 

X 


1 


(j) = 

F V ms 2 +G a s 2 +G v s+K 


(B38) 


These equations are made non- dimensio nal, as was done for the transmissibility functions, 
by div iding by the mass m of the isolated payload and defining the following: 0^= (K/m)^, £ = 

l/2Gy(l/Km)^, and Gg/m = b2/2nf 2 0 . Making the above substitutions gives, for velocity feedback. 


—( s ) = l 

F K v s 2 + 2£(D ) ,S + a>* 


(B39) 


for acceleration feedback. 



S 2 +- 


2 Vf c 


TS 2 +(o\ 


and for both acceleration and velocity feedback 


(B.40) 


-(s\ = r (B.41) 

F * 5* +^5-** +2^+40^ 

To show the effectiveness of the active systems, the ratio of X(s)/F(s) for an active system to 
X(s)/F(s) for die passive part of die system will be used. This ratio wifi be called die mobility 
effectiveness, Xj(s). Thus, if Xf is unity, the effectiveness of die active vibration isolation system in 
reducing force induced payload motion is zero, or die active portion of the system does nothing. If Xf 
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is zero, die effectiveness of die active potion of die system is complete and there is no motion of die 
isolated payload resulting from a finite applied force. If Xf is greater than unity, then die active 
portion of die vibration isolation system amplifies die effect of the applied force, giving rise to more 
payload motion with active feedback than without it. The equations for the effectiveness function for 
the different systems are, for velocity feedback, 

X (s) = £l±lilgLa£±J»l <B 42 ) 

f ’ ** + 2fyo m s + m* 

for acceleration feedback. 


v (A _il±25l®-£±®^ 

A / W -“ 


s 2 +-^rS 2 +(al 


(B.43) 


2vf: 

where = l/2c(l/Km)^, and for both acceleration and velocity feedback. 


*/M. = 


y 2 +2^,c),s+to^ 


s 2 +-r~brs 2 +2£fi> Jl j+co 2 


(B.44) 


2*f; 


In terms of frequency response, the vector length Xf(s) is |Xf(/H>)|. This value is obtained 

in die same manne r as was done for transmissibility. The equations for the effectiveness of the 
systems become, for velocity feedback. 



for acceleration feedback. 
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(B.46) 


and for bo* acceleration and velocity feedback, 





i- 

M 

1®J 
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i + -hJ 
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1/2 


(B.47) 


Figures B.5-B.7 represent the effectiveness of the force actuated vibration isolation systems for a 
passive damping ratio of 0.05, which is typical of spring elements made of steel. The figures show 
the effectiveness term |X / (yco)| as a function of the non-dimensional gain terms. 


B 3 Transmissibllity Evaluation of Inertial Feedforward Controller 


In order to demonstrate die advantages of isolating a payload by the proposed feedforward 
inertial controller design, a simple one degree erf freedom spring mass damper Systran, shown in 
Figure B.8, is solved were F s is a servo force proportional to die inertial position and velocity of the 
support structure, specifically F s is defined as; F t = ku f + Cli / . This force is achieved by referencing 

an ac tuato r to the appropriate integrals of an accelerometer, attached to die support structure of the 

mass. The equation of motion for this system is: 

/7tr+c(i-u)+£(x-u)+]^ =0 


(B.48) 






203 



Figure B.6: Inertial acceleration feedback mobility curves. 
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Figure B.7: Inertial acceleration and velocity feedback mobility curves. 



205 


and substituting the defined servo force into (B.48), 

mx+c(x-u+u')+k(x-u+u') = 0. (B.49) 

Using die following definitions, die viscous damping factor £ = where (Oq is the natural 

frequency of vibration for die system, to^ n = k/tn. Equation (B.49) can be put into vibration notation 
using die above definitions and the equation of motion becomes: 

x + 2^mJi+G) 2 u x = Kco^l 2 £g>„i^1 -4-J- (B30) 

If one assumes monochromatic motion of the base structure and that the inertial sensor has some 
phase and amplitude error, one can write the displacement, u, as Asin(cot), where A is the 
displacement amplitude. Therefore, die inertial sensor, not including constants of integration, would 
give an inertial displacement of die support structure equivalent to u' as A'sin(cot + $), where A' is the 
displacement amplitude detected by the accelerometer. If we take die Laplace transforms of x, u, and 
u\ including the ratios u'(t) / u(t ), and u'(t) / u(t), equation (B.50) becomes: 

X(s)(s 2 +2^to n s+©^) = U(s^to^l- 

(B.51) 

Taking the frequency response of the transfer function and calculating the magnitude of this function 
one arrives at the following equation, which gives the frequency response transfer function or 
transmissibility of die isolated payload to a harmonic base disturbance. 



(B.52) 
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Equation (B .52) depicts the feedforward inertial cancellation approach to offsetting fee attenuation 
profile of a second order spring mass damper system, hi the previous formulation the spring and 
damper elements were assumed constant However, the approach explored, and documented in this 
work, has the advantages of a closed loop active magnetic system, where the relative spring and 
viscous elements are variables dependent cm controller loop gains. The error terms in equation (B.52) 
model non-ideal input sensors and solves the transfer function giving die frequency shaped 
transmissibility curve shown in Figure B.9. These error terms can be viewed as shaping functions 
which could be introduced into file closed loop control as frequency weighted terms to give specific 
response performance dependent on excitation profiles or user requirements. Figure B.9 illustrates 
the advantage of inertial cancellation for a modeled Systran. The solid curve gives die frequency 
response of the single degree of freedom system with a viscous damping factor £ = 0.1, and a natural 
frequency of (Dq = 2n(0.6 Hz). The dashed line is comparing die same relative based system with die 
addition of inertial cancellation using an inertial s tiffne ss and damp in g amplitude ratio of 1.0 with no 
error assumed in the inertial position and a phase shift of 0.4 radians in the inertial velocity. 
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Figure B.9: Frequency response of inertial and relative based one DOF configurations. 




Appendix C 


One-Degree-of-Freedom Magnate Circuit Actuator Suspension 


The purpose erf this appendix is tt> support the summary of die attractive relative suspension 
system presented as background for the prototype fecdforward/Feedback isolation system. The 
formulation of the relative suspension equations for an isolation system are very s i milar to equivalent 
formulations for rotating magnetic bearing systems. Such derivations can be also found in many 
papers on the subject of magnetic bearings [e.g., 24]. 

In analyzing a one degree of freedom attractive relative magnetic actuator suspension 
system, as shown in Figure C.l, small motions are assumed about some crater position so a linearized 
approach can be taken. We also assume die flux levels in the core and armature of the 
electromagnetic circuit are below saturation. The total magnetic flux in an air gap [24] is 


Neglecting leakage and fringe effects in the magnet, die force exerted by die magne t is 

<b 2 

F = -^ — 

' 2 M 


(C.1) 


(C.2) 


Therefore, substituting the relation for the magnetic flux in an air gap gives the relation 

jV ^ 2 " 2 

p~ 


MNV 
2h 2 


(C.3) 


For both air gaps, die total force exerted is dins 

U AN 2 i 2 

F=2F ' = iT~ 

Then the ac tual gap h and current i, for die system shown in Figure C.1, have die form 

h = h 0 -lfi 


(C.4) 


(C.5) 
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Figure C.1: Physical description of relative one degree of freedom suspension. 
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i = i b + Ai (C.6) 

The negative sign occurs in h because as the mass of the system moves up, the clearance in the gap 
decreases. Therefore, substituting these relations for h and i into the force r el ati on for an 
electromagnetic circuit one arrives at die following equation: 


u an 2 (l+m) 2 

F= , V * (C.7) 

{K-ifl) 

By assuming both and i are small, equation (C.7) can be approximated by a binomial expansion 
about the bias current and the gap h. This then gives a linear relationship between the force F and 
i, making it possible to use a linear control scheme. The binomial expansion of equation (C.7) 
becomes 


H.ANHl ( llfi 

The static equilibrium equation for the one DOF system as shown in Figure C.l is 

mg^ = F<h 


(C.8) 


(C.9) 


Thus, 


F - mm8 2T 


l _2MJV^ 


hi 


Thai by setting i in the linearized force equation to zero, one arrives at 


F-m* 1 JKQ 1 - | 

F-mg--kfit,- h] + k . 


Therefore, the position stiffness is, from equation (C.ll), 


(C.10) 


(Cll) 


2\l 0 AN 2 i 2 b 

“ Li 


(CM) 
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A position dis pl ace me nt /j toward die magnet increases the force in dial same direction. 
An actual spring would apply a force tending to restore die initial position of die beam at Now 
setting in the l inea ri zed force equation to zero one gets a linearized force component proportional to 
die bias current. This has been called a current stiffness by many [24]; however, a change in i does 
not tend to restore die beam to its original position, but an increase in c u rr e nt will tend to force die 
beam away from its steady state value always towards the magnet The following relation 
demonstrates this proportionality constant which arises from linearizing die force around a bias 
current: 


„ 7 .... \l.AN‘i‘ h . 

F. mg --k l A t ^— j ^ + h , ■■ 

Therefore this proportionality constant becomes 

2M/n 

hi 

The dynamic equation for the system of interest can be written as 

T 2;2 


(C.13) 


(C.14) 


/ | i 0 AN 2 ii ... . dQ . 

mg — + AF = °— 2 b --+k i Ai+ml l -^T+ W 


(C.15) 


Since, from the static equilibrium case, F 0 was set equal to the force needed to support die beam at its 
equilibrium position. 


F 0 =mg— = 


l _ \i b anHI 


(C.16) 


2A hi 

then dynamic equation of this system around its static equilibrium position, which is set by the 
amount of bias current, can be represented by 

A F = kjAi+mli -^-+k 9 Ql 1 


(CAT) 



213 


By assuming the current i = + i c , where is the bias current and i c is the control current, then A i 

= i c . And assume that at equilibrium, or at 8=0, l^h^ The dynamic response equation becomes 

A F = kj e +m/,— +^0^ (C18) 

In order to stabilize die system one needs to take die derivative of die position signal and use 
both position and velocity to control the system. This portion of the control circuit can be m odeled as 
a position constant kg and a velocity constant kp These gains are adjustable and effect the s ti ffn es s 
and damping coefficients acting on die system. The transfer function of this portion of the circuit can 
be modeled as 



+ 


^2 S 

l + t 2 $ 


k t +(k t +k r )z 2 s 

1+t 2 s 


Then current amplification is described by a constant as follows: 


(C-19) 


— = k a (C20) 

A block diagram can be formed as shown in Figure C.2, using force as the input and position as 
output The transfer function becomes 


1 



m s 2 +k 9 

1+ - 1 

kikjk p \k t +(k g +* r )T 2 s] 

ms 2 +£3 

(l+t 2 s)(l+x,j) 

b J 


(C21) 


In order to determine die theoretical stiffness and damping coefficients, die dynamic 
equations of the one DOF system are written and compared to a lumped second-order model. By 
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Figure C2: Control block diagram. 
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using the Wock diagram. Figure C.2, an equation can be written involving position 0 and control 
current i c as follows: 


. )v] 

,e (i+vXi+v) 

By using the relation 

AF = k i i c +(ms 2 +k 9 )e 

this system of equations can be put into matrix form as follows: 


(C22) 


(C23) 


TAFT 

ms 2 

T = 

k a k p [k t +(k g +k r )x 2 s] 

L u J 

(I+VXI+V) 



(C.24) 


Then in solving this matrix formulation of the system dynamic equations for 6, one arrives 
at die following equation for 6(s): 


where 


Thus 


e(s)= 


AF 


(m+K)+ (llvXl+v) 1 


AF = k,-i c + (ms^ + k q)6. 


(C.25) 


(C26) 



k,i e +(ms 2 +k t )Q 


(i+vXi+v) 


(CJ27) 
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Setting this relation equal to the lumped second-order system of a spring dashpot configu ration one 
can get the theoretical stiffness and damping values for die closed-loop magnetic circuit 
configuration: 


©0g>) = - 


+(*.-!« o 2 )0 


-mcol+ 


)■ 


k i k.k f [k t +{k t +k r )xJ(o 


{l+TjaMl+xJa) (C28) 

/(o>) 

(k tq -nm 2 )+ jdK^ 

Thus, from the equality given in equation (C.28), one obtains die following r elatio ns: 

j {<& ) = kji' c + (kg -mw 2 ^, (C.29) 


and 


= ite+^ 


*,*.*,[*« +(*, +k r )x 2 j(i>] 

{l+x 2 j(o)(l+xJ(o) 


(C.30) 


Jtrv 

f — < 

W, 

*, +(* f +/: r )x 2 yti>]] 

C €0 

* CO 

(1+t 2 >0))(1+Vcd) j 


(C.31) 


Therefore, in evaluating the equivalent stiffness and damping for die closed-loop system one 
comes up with the following relations: 


wp 2 )+(*, +*,X*2 +*i 
49 + (l-TjT^J+Cx 2 +i,)V 

C * (I-TjTjCo’J+OCj+tJV 



(C.32) 


(C33) 



Appendix D 


Transfer Function Determination 


This appendix will describe how th e signal analyzer used extensively in the design and 
development erf the active isolation systems determines the transfer function of signal c ond i tion i ng 
filters or any other system. The generation of random noise by the signal analyzer is described first 
and then the processing accomplished by the analyzer is discussed. 

The si gnal analyzer generates an exponentially time-correlated random process with a 
correlation time, T, and mean-square value of o^. This process is effectively Gauss-Markov, which 
imp lies wide-sense stationarity and ergodicity, and therefore has properties which make the 
mathematics tractable. The autocorrelation function, RjmCt), of tins process, n(t), is shown in Figure 
D.l and can be described by 

R an M = oVVl'l. 



Figure D.l: Autocorrelation of an Exponentially Time-Correlated Process 
The analyzer obtains the antospectral (or power spectral) density of this random process ty 
takin g the fast Fourier transform of the autocorrelation function. In the frequency domain, the 
antospectral density, G^©), of the random process, n(t), is illustrated in Figure V>2 and is given by 
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G„((0) = 


2a 2 /r 

(q 2 +(1/T) 2 ' 


(D.l) 



Figure D.2: PSD of Exponentially Time-Correlaled Process. 

The result is a band-limited process, displaying equal power at all frequencies within die band of 
interest The half-power point occurs at 1/T, which corresponds with die highest frequency selected 
through the analyzer span control. Through judicious selection erf the frequency span, a random 
process displaying die desired characteristics is obtained. 

For analytical purposes, white noise as opposed to wideband noise, is desirable to work with 
for a number of reasons. It is easily manipulated with engineering tools and the theory is complete 
and practical. Whiteness implies that the noise value is not correlated in time and dial the noise has 
equal power in all frequencies. Since this would result in a noise with infinite power, a white noise 
process does not realistically exist Although this concept does not exist it is useful to consider. Any 
physical system has a certain frequency bandpass, Le., a frequency range of inputs to which it can 
respond. Above this range, the input either has no effect or the system attenuates the effect and 
essentially the input does not exist beyond a certain stopband. 

Typically a system will be driven by wideband noise, one having constant power at 
frequencies within the passband. On this plot, a white noise would merely extend this co nstant power 
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level across all frequencies. Within the bandpass of the system, the fictitious white noise looks 
identical to wideband noise. 

In some instances, the noise power level is not constant over all frequencies within the 
system bandpass, or die noise may be time correlated. For such cases, a white noise put through a 
small linear filter can duplicate virtually any form of time-correlated noise. This is the case of die 
random process generator in the signal analyzer. 

To obtain a random process exhibiting die desired autocorrelation and antospectral density, 
the analyzer drives a first order lag system with a zero mean, white Gaussian noise of strength Q, 
where Q = 2o^/T. After applying the white noise to the lag circuit and allowing transients to die out, 
a steady state process with die desired statistics results. 

It may be obvious why it is necessary for the analyzer to alter the half power point of the 
random process it generates, when alternatively it could generate a single process with an 
exceptionally wide band. In actuality, the analyzer can only generate a process that makes die 
integral of the PSD constant. Reducing the half power point (increasing the correlation time) results 
in a higher power level distributed over die frequency band. likewise, increasing the half power 
point, decreases the power level. When exciting a system across its bandpass, it is desirable to 
distribute the highest amount of power within die frequency band of interest Therefore, it is 
necessary for the analyzer to provide for die adjustment of the random process statistics through the 
selection of die proper frequency span value. 

The data processing performed by the analyzer in seeking to determine a test item's transfer 
function requires the basic setup shown in Figure 3.13. The random noise process is injected into the 
system and sampled by a channel of the analyzer. The system output is sampled by a second channel. 
Prior to digitization, both channels of data are filtered to prevent aliasing. Because processing is 
accomplished in real time, biases and trends are not removed. Using the data collected from the 
random process, x(t), the autocorrelation function Rxx(t) is calculated. The random process data, 



220 


together with the data collected from the system output, y(t), allow the calculation of the cross 
correlation function Rxy(r). Both of the correlation functions are then fast Fourier transformed, using 
side love suppression techniques, to give die autospectral density, Gxx(to), and the cross spectral 
density, Gxy(co). Using complex forms of these densities, die system transfer function is computed by 


ff((D) = 


Gjca) 


Oja)' 

The magnitude and phase result is displayed on the analyzer screen. 


(D.l) 


D.l Autocorrelation and Crosscorrelation Functions 


Given an accelerometer output measuring the random vibrations of some structure, the 
voltage output from the accelerometer, x(t) will be a random signal. If two voltage measurements, vj 
and V2 are taken at two different times, tj and t j + 1 respectively, and this process is repeated many 
times giving multiple pairs of measurements {vj^}, where in each pair, vj is measured i seconds 
before V2- If the time interval t is extremely small, then vi and V2 are statistically likely to be nearly 
the same value. 

For most random processes, an increase in die time duration ? results in a lower probability 
that vj and V2 became less interdependent, or correlated as die time interval ?, is increased. The 
autocorrelation function may be described as a measure of die correlation between two successive 
measurements, as a function of the time elapsed between these two measurement. It follows that the 
autocorrelation function is also an indication of how rapidly the function is changing with time. If 
the function changes very rapidly, then the time duration t, will have to be extremely small in order 
to obtain good correlation between the voltage measurements vj and V2- 

The expressions for the autocorrelations of two random processes x(t) and y(t) are shown 
below for the special case in which the process means p^ and Py, and the autocorrelation functions 
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themselves are not changing with time. If this is not die case, die autocorrelation functions beco me 
functions of both time t, and the time inter val between measurements, t. It is assumed that the 
process means and die autocorrelation functions process for a certain length of time do not change. 
Such processes are often called stationary processes. It is also assumed in the expressions below that 
the random processes are both zero mean. 

RM = E{x{t)x(t + t)} = + i)dt (D.2) 

Ry,W = B{y(t)y(t+x)} = fyfiy{t+%)dt (D.3) 

In these expressions die aerator E{ } is used to denote die arithmetic mean of the set of values 
mr.i nsftd by the brackets. For example, in the expression for R^fr), all elements of die time history 
x(t) which form pairs separated by a fixed time interval i, are id e nt i fie d. The elements of each pair 
are multiplied together, and the arithmetic mean of all of these products is then calc ulat ed. The 
process is repeated for the next value of t, until all values of t have been exhausted. 

The crosscorrelation function is closely related to die autocorrelation function. The only 
difference between the two functions is die source of the measurement V 2 in the above discussion of 
die autocorrelation function. The first voltage measurement vj is a measurement of die random 
process x(t) as before, but die second measurement vj, is a measurement of some other random 
process. y(t). As before, die two measurements are taken within t seconds of one another. The 
expression for the crosscorrelation function between x(t) and y(t) is shown below. As was assumed in 
die presentation of the autocorrelation functions, the two random processes are zero mean. In 
addition, the process means and die crosscorrelation function are not changing with time. 

Rjr) = E{x(t)y(t+%)} = (D.4) 
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D.2 Autospectra and Aotospectral Density Functions 

A random process which has a mean and an amocondation function that are not changing 
with time is said to be stationary. The signal content in die frequency domain, of a stationary random 
process is given by die Fourier transform of the autocorrelation function as 

Sjjl o) = ~ J “ CD S) 

where: j = imaginary operator = V— I, and co = frequency in rad/s. S^Qei) is called die two sided 
autospectrum of the random process x(t). This function is said to be two sided because it is symmetric 
about the © = 0 axis. The spectrum equation is modified by dividing by the bandwidth of interest, 
and by changing the limits of integration to reflect this bandwidth. If die signal has units of volts, 
then the division by the bandwidth gives die spectrum units of V^/(iad/s). The modified spectrum 
having the units of V^/(rad/s) is a measure of the density of signal per unit frequency increment. For 
this reason, the modified antospectxum is called an autospectral density function. The limits of 
integration in an autospectral density function must also be consistent with the bandwidth. For 
example, if an autospectral density is to be calculated for a signal over the frequency band ©^ 5 © S © 
f, then the limits of integration are from coj to ©f . Notice that the resulting function is not symmetric 
about the © = 0 axis. This function, defined over the above frequency band is said to be erne sided. In 
order to reflect all the energy in the signal over to die right hand side of the © = 0 axis, die spectrum 
must also be multiplied by 2. It should be noted dial die autocorrelation function R^Ct) is also a 
symmetric function. When die band limited aotospectral density function is multiplied by 2, only die 
portion of Rjocft) corresponding to positive values of x need be included. The complete, one sided 
autospectral density function, which is the standard for describing spectra is defined below. 


where bw = bandwidth (rad/s). 
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D 3 Crossspectra and Crossspectral Density Functions 

The one sided crosspectral density function between the random processes x(t) and y(t) is 
given below. This function differs from die one sided autospectral density function in that the 
crosscorrelation function Rjy(t) is used radio’ dim die autocorrelation function R n (t). 



Appendix E 


Input and Output Spectra for Multi-Dimensional Transfer Function 
Calculations 

The following figures are given as reference to the transfer functions calculations from 
Chapter 4 for the multi-dimensional Learjet transfer functions ealmlan»H for the relative and inertial 
control cases. Figures E. 1 through E.4 (a) and (b), give the input and output spectral densities for the 
X direction for the single input reference locations one through four for both trajectories, relative and 
inertial, respectively. Figures EL5 through E.8 (a) and (b), give the input and output spectral density 
functions for the Y direction, and Figures E.9 through E.12 (a) and (b) give the input and output 
spectral density curves for the Z direction, for the single input reference locations one through four 
for both relative and inertial control cases, respectively. Figures E.13 and E.14 (a) and (b) give the 
input and output spectra for the 0 X rotation for the trunnion location and the aircraft location, 
respectively. Figures E.15 and E.16 (a) and (b) give the input and output spectra for the rotation 
as defined in equations (4.17) and (4.18), and Figures E.17 and E.18 (a) and (b) are the input and 
output spectra for the rotation as defined in equations (4.19) and (4.20). These spectral density 
functions have been plotted with the same minimum and maximum values for both relative and 
inertial control cases. 
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(a) Payload under relative control. 

Figure E.1: Output and input power spectra for platform x and aircraft location 1 x. 
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(b) Payload under inertial control. 

Figure EJ.: Output and input power spectra for platform z and aircraft location 1 x. 
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(a) Payload under relative control. 

Figure E 2 : Output and input power spectra for platform x and aircraft location 2 x. 
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(b) Payload under inertial control. 

Figure E 2 : Output and input power spectra for platform x and aircraft location 2 x. 
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(a) Payload under relative controL 

Figure EJ: Output and input power spectra for platform * and aircraft location 3 x. 
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(b) Payload under inertial controL 

Figure EJ: Output and input power spectra for platform x and aircraft location 3 x. 
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(a) Payload under relative controL 

Figure E.4: Output and input power spectra for platform x and aircraft location 4 x. 
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(b) Payload under inertial control. 

Figure E.4: Output and input power spectra for platform x and aircraft location 4 x. 
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(a) Payload under relative control. 

Figure ELS: Output and input power spectra for platform y and aircraft location 1 y. 
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(b) Payload under inertial control. 

Figure EJ : Output and input power spectra for platform y and aircraft location 1 y, 
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(a) Payload under relative control. 

Figure E.6: Output and input power spectra for platform y and aircraft location 2 y. 
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(b) Payload under inertial control. 

Figure E.6: Output and input power spectra for platform y and aircraft location 2 y. 


.1000 1.000 10.00 7S.00 

Log Frequency (Hz) 



Hu to spectrum 


237 



(a) Payload under relative control. 

Figure E.7: Output and input power spectra for platform y and aircraft location 3 y. 
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(b) Payload under inertial control. 

Figure E.7: Output and input power spectra for platform y and aircraft location 3 y. 
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(a) Payload under relative control. 

Figure E.8: Output and input power spectra for platform y and aircraft location 4 y. 
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(b) Payload under inertial control. 

Figure EJ5: Output and input power spectra for platform y and aircraft location 4 y. 
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(a) Payload under relative control. 

Figure IL9: Output and input power spectra for platform z and aircraft location 1 z. 
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(b) Payload under inertial control. 

Figure EL9: Output and input power spectra for platform z and aircraft location 1 z. 
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(a) Payload under relative control. 

Figure E.10: Output and input power spectra for platform z and aircraft location 2 z. 
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(b) Payload under inertial control. 

Figure E.10: Output and input power spectra for platform z and aircraft location 2 z. 
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(a) Payload under relative control. 

Figure E.11: Output and Input power spectra for platform z and aircraft location 3 z. 
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(b) Payload under inertial control 

FlgureE.il: Output and input power spectra for platform z and aircraft location 3 z. 
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(a) Payload under relative controL 

Figure E.I2: Output and input power spectra for platform z and aircraft location 4 z. 
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(b) Payload under inertial control. 

Figure E.12: Output and input power spectra for platform z and aircraft location 4 z. 
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(a) Payload under relative control 

Figure E.13: Output and input power spectra for platform G x and trunnion 0 r 
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(b) Payload under inertial control. 

Figure E.13: Output and input power spectra for platform (^and trunnion 0 r 
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(a) Payload under relative control 

Figure E.14: Output and input power spectra for platform 0 x and aircraft 9 r 
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(b) Payload under inertial control. 

Figure E.14: Output and input power spectra for platform 0 X and aircraft 0 r 
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(a) Payload under relative control. 

Figure E.15: Output and input power spectra for platform and aircraft ^ 2 - 2 - 
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(b) Payload under inertial control. 

Figure E.15: Output and input power spectra for platform and aircraft 4y2-3< 
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(a) Payload under relative control. 

Figure E.16: Output and input power spectra for platform ^ and aircraft 
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(b) Payload under inertial controL 

Figure E.16: Output and input power spectra for platform and aircraft 
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(a) Payload under relative control. 

Figure EJ7: Output and input power spectra for platform P z and trunnion P z . 
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(b) Payload under inertial control. 

Figure EJ.7: Output and input power spectra for platform P z and trunnion P z . 
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(a) Payload under relative control. 

Figure E.18: Output and input power spectra for platform p z 211(1 aircraft p z « 
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(b) Payload under inertial control. 

Figure E.18: Output and input power spectra for platform f) z and aircraft f) z . 
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